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Abstract 



The aim of this thesis is to derive new gradient estimates for parabolic equations. The 
gradient estimates found are independent of the regularity of the initial data. This 
allows us to prove the existence of solutions to problems that have non-smooth, con- 
tinuous initial data. We include existence proofs for problems with both Neumann and 
Dirichlet boundary data. 

The class of equations studied is modelled on mean curvature flow for graphs. It 
includes anisotropic mean curvature flow, and other operators that have no uniform 
non-degeneracy bound. 

We arrive at similar estimates by three different paths: a 'double coordinate' ap- 
proach, an approach examining the intersections of a solution and a given barrier, and 
a classical geometrical approach. 
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Chapter 1 

Preface 



Mean curvature flow 

Let {M t } be a family of hypersurfaces, each smoothly embedded in R n+1 and indexed 
by t e [0,T]. We say M t is moving by mean curvature flow when 



where H is the mean curvature vector at x e M t . 

In the last twenty-five years, this flow has been the subject of concerted study, as 
have other geometric flows such as the Ricci flow and Gauss curvature flow. Notable 
results have included Grayson's proof that the curve shortening flow (which is mean 
curvature flow, reduced to one space dimension) shrinks embedded curves to a spher- 
ical point [1 7] and Huisken's proof that convex surfaces become spherical under mean 
curvature flow (18) . 

If we observe that H = A Mt x, where A Mt is the Laplace-Beltrami operator on the 
manifold M t , then it seems natural to consider mean curvature flow as the heat flow 
for manifolds. 

Unlike classical heat flow, this is a nonlinear operator, but it still has some of the 
same attributes; in particular, this flow exhibits a smoothing property. In fj4|, Ecker 
and Huisken showed that if the initial surface is given by a locally Lipschitz graph, then 
there exists a smooth solution for positive times. In this thesis, this result is extended 
to non-Lipschitz initial conditions. 



Mean curvature flow and parabolic differential equations 

Chapter |2] is a short introduction to parabolic differential operators and some key re- 
sults in mean curvature flow. 

If M t is locally represented as graphs for some u : 9, x [o,T] — ► R, then u will 
satisfy the parabolic equation 



du 
~dt 



y/l + \Du\ 2 div 



y/l + \Du\ 2 



Du 



) 




1 + \Du\ 2 



DiuDju 



) 



Dau. 
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This places mean curvature flow in the setting of classical parabolic partial differential 
equations, the framework for most of this thesis. 

Motivated by this setting, we examine other parabolic operators that have similar 
diffusion properties to mean curvature flow. This includes anisotropic mean curvature 
flow, a generalization of mean curvature flow arising in many physical applications. 



Gradient estimates using a 'double coordinates' approach 

In this thesis three distinct methods are used for deriving gradient estimates. The first 
of these is introduced in Chapter|3l It originated with Kruzkov in |22). Given a parabolic 
equation in one space dimension 

u t = F (u xx ,u x ,u,x,t) 

one can form an evolution equation for the difference w(x, y, t) = u(x, t) - u(y, t), 

w t = F (u xx ,u x ,u(x),x, t) - F (u yy ,Uy,u(y),y,t) . 

Under favourable conditions on F, one can use the maximum principle and an appro- 
priate barrier to find estimates for w that depend on \x - y\. Letting \x - y\ — > will 
then give a gradient estimate for u. 

In this thesis, Kruzkov's method is extended by making use of the full Hessian 

[D 2 w] 

rather that just the diagonal elements w xx and w yy . In the one-dimensional case this 
is of little importance, but in the higher-dimensional case it gives us greater scope to 
choose barriers. 

In Chapter|3l we also describe a barrier which begins with an unbounded gradient, 
but instantly becomes smooth. 

Such barriers allow estimates that are independent of initial gradient bounds, and 
that therefore may be used to prove the existence of solutions to parabolic equations 
with continuous initial data. This is one of the main features of the gradient estimates 
in this thesis. 

The 'double coordinate' method is extended to higher dimensions in Chapter|5]for 
a class of operators that have similar diffusion to the mean curvature flow. 

As this class includes anisotropic curvature flow , this is a significant improvement 
to the existing regularity theory. 

Gradient estimates are found for both entire periodic solutions and boundary value 
problems. 



Existence results 

The gradient estimates of Chapters|3]and|S]may be used to extend standard existence 
results. 



Chapter 1: Preface 
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We do this for a class of parabolic equations in the one dimensional case in 
ChapterEl for mean curvature flow with Neumann boundary conditions in Chapter 
IH and for mean curvature flow with Dirichlet boundary conditions in Chapter|7| 

Gradient estimates by counting intersections 

The second technique for finding gradient estimates is found in Chapter|8l and it in- 
volves examination of the intersections between a given solution u and a barrier ip. 

In 0, Angenent proved that the number of points in the zero set — the set where 
u(x,t) = — of a solution to a parabolic equation in one space dimension is non- 
increasing. 

This is applied to the difference u - p. The intersections of u and <p are the zeroes 
of u - ip, and so Angenent's results allow us to show that u and p do not develop new 
intersections as they evolve. 

When two functions intersect only once, the gradient of one of them dominates the 
gradient of the other at that point. Tailoring the barriers gives us estimates for u x (x, t) 
in terms of the height of u(x,t), the time t, and (in the case of bounded domains) the 
distance of x from the boundary. Again, there is no dependence on an initial gradient 
estimate. 

As Angenent's results are limited to equations in one space dimension — it is 
difficult to imagine what a generalization of these results to higher dimensions would 
look like — this technique applies only to parabolic operators in one space dimension. 

These methods will apply to a wide range of parabolic operators, provided that 
suitable barriers exist. In particular, we can apply this method to the class of operators 
studied in Chapter^] 

Gradient estimates using a geometric approach 

The third method for finding gradient estimates, in Chapter|U is a rather geometrical 
approach found in the classic Ecker-Huisken curvature flow papers [T3], [T4|, [18], 
ITSl and (20). 

A maximum principle is applied to the difference Z = v - <p, where v is a "gradient 
function" (for example, v = y/1 + \Du\ 2 ), and p is a barrier. In contrast to the earlier 
two approaches, this creates a direct estimate for the gradient Du itself, rather than for 
the difference u(x) - u(y). 

We apply this to the mean curvature flow (re-creating some of the results obtained 
in earlier chapters) and also to the anisotropic mean curvature flow, under some re- 
strictions on the degree of anisotropy allowed. Results for entire periodic solutions 
and strictly interior results are found in both cases. The estimates found are again 
independent of initial gradient bounds, but dependent on the height. 

Appendices 

An inventory of standard results for parabolic partial differential equations, relevant to 
the existence results in Chapters 01 16] and is included for the convenience of the 
reader. There is also a nomenclature listing. 



Chapter 2 



Mean curvature flow and parabolic 
equations 



matrix of derivatives F~ 



2.1 Parabolic equations 

An operator P:Ix § n xl n xlxO^lx[0,T]is considered parabolic on a domain 
when 

P(z + a,r + r),p,q,x,t) > P(z,r,p,q,x,t) 

for any positive definite 77 e § n , positive number a, and any (z,r,p,q,x,t) in the do- 
main. Here, S n is the set ofnxn symmetric matrices. 

In this thesis we look only at operators of the form 

Pu = P(—ut, D 2 u, Du, u, x, t) = —ut + F(D 2 u, Du, u, x, t). 

If F is differentiable with respect to the first variable, then P will be parabolic if the 

" dF 1 

- — is positive definite. 

If we can write F(r,p,q,x,t) = a ij (p,q,x,t)rij + b(p,q,x,t) for some symmetric 
a:M n xlxfix[0,T]^ then we call the operator quasilinear. It is to quasilinear 
operators that we will pay most attention in the following pages. 

In this case, F r = [a lj ] and so the operator is parabolic on S if [a ij q, x, t)] is 
positive definite for all (p, q, x, t) e S, where S is a subset of 1" x K x O x [0, T]. 

We denote the maximum and minimum eigenvalues of [a ij (p, q, x,t)] (or F r ) by 
\(p, q, x, t) and A(p, q, x, t). 

If the ratio A/A is bounded on S, then P is called uniformly parabolic on S. 
An operator is parabolic with respect to a function u when P(D 2 u,Du,u,x,t) is 
parabolic. 

When A ^ — for example, when a ij (Du, u, x, t) = \Du\ p ~ 1 5 ij , as in the parabolic 
p-Laplacian equation — such an operator is called degenerate. 

The key to much of the theory used here is the Comparison Principle. As presented 
in |25l: 
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Theorem 2.1 (Quasilinear comparison principle 1). Suppose thatP is a quasilinear 
parabolic operator 



Let u and v be in C 21 (Q \ VQ) n C(Q) and let P be parabolic with respect to either u 
orv. Then if Pu > Pv in Ti \ Vfl and ifu<v on V, then u <v in 0. 

Here V denotes the parabolic boundary 



The proof of this theorem is simple, and is an excellent illustration of later arguments. 

Proof: Suppose that there is an interior point x at time t > where u = v for 
the first time. Since this is an internal maximum of w = u — v, Dw = Du(x ,t ) - 
Dv(x ,t ) = and Dijw = Diju - Dijv must be negative semi-definite. Now, 



However, as this is the first such maximum we must have w t > 0, which give us a 
contradiction. It follows that u < v. □ 

We also use the following form, again as in [25J : 

Theorem 2.2 (Quasilinear comparison principle 2). Suppose thatP is a quasilinear 
parabolic operator 



and that there is an increasing function k(M) such that b(p, q, x, t) + k(M)q is a de- 
creasing function ofqonQx [-M,M] x M. n for any M > 0. Ifu andv are functions in 
C 2 ' 1 ^) n C(0) with Pu > Pv in Tl\VQ. and u < v on vn, and if P is parabolic with 
respect to either u orv, then u<v in Tl. 



2.2 Mean curvature flow 

If our family of hypersurfaces M t is also a family of embeddings F t -. M n — > M. n+1 , then 
we can write mean curvature flow as 



Pu = —ut + a 13 (Du, u, x, t)DijU + b(Du, u, x, t). 



r(nx [o,r]) := n x {o}udn x [o,r]. 



w t = u t - v t 



= —Pu + a l] (Du, u, x, t)DijU + b(Du, u, x, t) 

+ Pv — a lJ (Dv, v, x, t)DijV — b(Dv, v, x, t) 
< a 13 (Du,u,x,t)DijW. 



Pu = —u t + a 13 (Du, x, t)DijU + b(Du, u, x, t), 




(2.1) 



where H is the mean curvature vector at F t (x) e M t . In the case that M t can be 
written locally as a graph over a set Vt e W 1 , we write F t (x) = (x,u(x,t)), and can 
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calculate geometric quantities such as the upwards unit normal 

(-Du,l) 



the metric on the surface 
the second fundamental form 



" + \Duf 

9ij = $ij + DiiiDjU, 



hij 



DijU 



71 + \Duf 



and the mean curvature 

H = g lj hij = - div 



Du 



The mean curvature vector is H = Hu, and so (if we remove movement tangential 
to the surface) mean curvature flow for graphs is given by 

du r 2 / Du \ ( DiuDju \ 

— = V 1 + \Du\ 2 div - = 5ij - \ Diju. 

dt \y/l + \Du\ 2 J \ l + \Du\ 2 J 

In the case when n = 1, this reduces to curve-shortening flow 

dll U Tr 



dt l + u 2 x 

With reference to the previous section, note that the largest and smallest eigenval- 
ues for mean curvature flow for graphs are A = 1 and A = (1 + H 2 ) -1 , so it will be 
uniformly parabolic only when the gradient is bounded. 

Whether studied in a geometric setting as \2A) . or as a special case of a quasilin- 
ear parabolic differential equation as i|2.2| . the comparison principle has been crucial. 
Applied to mean curvature flow, it gives: 

Theorem 2.3. Let M t and M[ be two smooth compact surfaces moving under mean 
curvature flow. If they are disjoint at the initial time, they are disjoint at later times. 

We can also make similar comparisons between surface s with bou ndaries, and be- 
tween other quantities (such as the gradient function v = y/1 + \Du\ 2 ). The following 
theorem from (TU is one such result. 

Theorem 2.4 (Interior gradient estimate). Suppose thatu satisfies the mean curva- 
ture flow equation (|2.2| on a cylinder B R (y ) x [0, T]. Then we have an estimate for 
the gradient at the center of the ball at later times: 



y/l + \Du(y ,t)\ 2 

, h 4- 1 Dnin nTis pyn 

R 2 



< Ci sup yjl + |£>u(y,0)| 2 exp 



/ \ 2 

C 2 I \ 

sup u(y,t) - u(y ,t) 

Bfl(j/ )X[0,T] J 



Chapter 2: Mean curvature flow and parabolic equations 
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where C\ and C 2 depend only on n. 

We will also use the following a priori estimates for higher derivatives, from the 
same paper: 

Theorem 2.5 (C 2 interior estimate for mean curvature flow). Suppose thatu satis- 
fies (ED on B R ( yo ) x [0, T]. Then for arbitrary < 6 < 1 the estimate 

sup \A\\t)<c{n){l-e 2 r 2 (^ + \\ sup {l + \Du\ 2 ) 2 

Bg R (yo) \ U V B R (y )x[0,t] 

holds for alio <t<T. Here \A\ 2 = h tj h kl g lk g kl . 

Theorem 2.6 (C k interior estimate for mean curvature flow). Suppose that u sat- 
isfies \2.2) on B R (y ) x [0,T]. Then form > and arbitrary < 9 < 1 we have the 
estimate 

sup \V m A\ 2 (t)<c m (l-9 2 r 2 U + - 
B eR (yo) \ K 1 

where c m is a constant depending on n, m and sup^^^ [ 4 ] (1 + iDul 2 ) 1 / 2 . 

A bound on \A\ gives a bound on \u\ C 2, since (using coordinates in which h is 
diagonal) 

|^| 2 = h ij h kl g ik g kl 
= h u g ik g kl h u 



n„m2 J2( hi 



- {l + \Du\ 2 ) 2 

- (TT^2)3 s S p|A ' n|2 ' 

where we have used that the smallest eigenvalue of g ij is (1 + \Du\ 2 )~ l . So, 

\D ijU \ < {l + \Du\ 2 f' 2 \Al 

and in a similar manner, bounds on derivatives of \ A\ give bounds on higher derivatives 
of u. 

These estimates may be used to show long-time existence results, such as the 
following from [13] : 



Theorem 2.7. If u is a locally Lipschitz, entire graph overW 1 , then there is smooth 
solution to E3 for all t > 0. 



In the following pages, we will derive new existence results of this sort. 



Chapter 3 



Gradient estimates for parabolic 
equations of curve shortening 
flow type in one space dimension 

In this chapter we outline gradient estimates for a class of parabolic equations in one 
space dimension. 

This chapter takes inspiration from the work of Huisken in |2TJ, where he inves- 
tigated embedded plane curves evolving by curve shortening flow by looking at the 
evolution equation for the quotient of d(p, q), the distance between two points p and q 
in the metric of the plane, and l(p, q), the length of curve between p and q. This intro- 
duced a double set of space coordinates (those around the point p and those around 
the point q). At a maximum point, the first and second derivative conditions give strong 
conditions at both p and q, allowing close examination of all possible situations. An 
application of the maximum principle resulted in a new proof of Grayson's theorem 
regarding the evolution of embedded curves. 

In this chapter, we follow the approach of Kruzkov in (22) (and well described in 
Lieberman's book [25J, chapter XI, section 6). 

If u(x,t) solves a parabolic partial differential equation in one space variable, then 
v(x,y,t) = u(x,t)-u(y,t) solves a parabolic equation in two space variables, for which 
we can seek a barrier. 

In the paper cited, Kruzkov was interested in fully nonlinear equations 

u t = F(u xx ,u x ,u,x,t) 
with uniform parabolicity condition 

d 

—F(r,p,q,x,t) > A > 0. 

In this section, we do not require uniform parabolicity, but in order to show existence 
of the barriers, we will require a scaling similar to that of the curve shortening flow 
equation 
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in that — ~ \p\~ 2 for large \p\. 

We begin with a description of the ideas motivating the method. The notation ' 
will indicate derivatives with respect to the space variable, which I hope I will use only 
where this is unambiguous. 



3.1 Outline of the 'double coordinate' method 

Consider a smooth h:Mx[0,T)^K satisfying 

u t = a(u x , u, x, i)u xx + b(u x ) 
u(x, 0) = Uq. 

Let Z:Mxlx (0, T) -> R be given by 

Z(x,y,t) = u(y,t) - u(x,t) - <p{\y - x\,t), 

where <p is some smooth function. 

Suppose now that Z attains a maximum at some point (x,y,t), with y > x. At the 
maximum point, the first derivatives are zero, and so 



= Z x = -u'(x, t) + <j)'(y - x, t) 
= Z y = u'(y, t) - <p'(y - x, t). 



(3.1) 



Similarly, the matrix of second order partial derivatives is non-positive, by which we 
mean that for all v e R 2 , v T [D 2 Z]v < 0, where [D 2 Z] is the Hessian matrix 



^11 


Z x y 




Zyx 


Zyy_ 





-u"(x,t) - <j)"(y - x,t) 4>"(y-x,t) 

(j)"(y-x,t) u"{y,t)-(j)"(y-x,t)^ 



If we now consider the evolution equation satisfied by Z, 

dZ 

— = u t (y,t) - u t (x,t) - 4> t {\y - x\,t) 

= a (u y , u(y, t),y, t) u yy (y, t) - a (u x , u(x, t), x, t) u xx (x, t) 
+b{u y ) - b(u x ) - 4> t (\y - x\,t); 

and if we take this at the local maximum we have 

dZ 



(3.2) 



dt 



trace 



i>\ u (y, t), y, t) (Zyy + 0") - a (</>', u ( x , t)i x i t) {-Zxx - <P") 

a((f>',u(x,t),x,t) ci 

c 2 a(<t>' \u(y,t),y,t)_ 

+ a ((/)', u(y,t),y,t) <f>" + a (</>', u(x, t), x, t) <f>" - (ci + c 2 )4>" - 4>t, 





Z X x 


Z X y 




Zyx 


Zyy_ 



for some c\ and c 2 . If the first matrix above is positive semi-definite, then as [D 2 Z] is 
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negative semi-definite, the trace above is non-positive and 



8Z_ 

dt 



< 4>" [a ((/>', u(x,t),x,t) + a ((j/,u(y,t),y,t) -c\- c 2 ] - <f> t . 



A useful choice for c\ and c 2 that makes the first matrix positive semi-definite is c\ = 

-a(<j)',u(x,t),x,t), c 2 = -a((p',u(y,t),y,t); then 



The idea now is to choose <fi in a way so that at the local maximum, Z t < 0. We 
begin by observing that for simple equations, a solution to a simplified version of the 
equation itself is acceptable for use as the barrier <f>. 

Remark: We could simplify the method by choosing c\ = c 2 = 0, in which case the 
factor of 2 is absent from Q3.3fr . The use of cross-derivatives will be important when 
we extend this method to higher dimensions. 



3.2 An estimate for periodic solutions 

Theorem 3.1 . Suppose u : R x [0, T) -> R is a H 2 , periodic and bounded solution of 



8Z_ 

dt 



< 2 (a (<//, u(x,t),x,t) + a ((/)', u(y,t),y,t) )(f>" - fc. 



(3.3) 



ut = a(u x )u xx + b(u x ) 



(3.4) 



with initial condition 



u(-,0) = u 



withu(x,t) = u{x + L,t) and osc u(-, t) < M. 
If ip : R x (0, T) -> R is a solution of 



ft > a(f x )f; 



XX ) 



with the initial and boundary conditions 



ip(x, t) -> 1 as t -> for x > 0, 

9?(0,t) = fort > 0, 

V9(x, t) — > l as x — >■ oo for i > 0, 



fr/en 




Proof: Following on from the previous remarks, we set 



Z(x,y,t) := u(y,t) 



u(x,t) - 4>(\y - x\,t) 



and choose <j>(z,t) = M(p(z/M,At/M 2 ). 

As t 0, ^ M for all 2 / 0, and so Z(x, y, 0) < 0. 
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As u is periodic, Z(x, y, t) = Z(x + L,y + L, t) and so Z is periodic over strips 

{ (x, y) : 2nL < y + x < 2(n + 1)L } . 



Within each strip, Z(y,y,t) = and 

Z(x, y, t) = u(y, t) - u(x, t)-(p<M-<j)^0 

as \y-x\ — > oo, so Z attains its spatial maximum in 
each strip, and hence in the entire plane, for each 
t > 0. We can calculate 

A - 4 ( 2 M 

> ^ w 

= HW, 

and in particular at a maximum point (x,y,t) with 
x 7^ y and Z non-negative, equation \3.3) becomes 

























\jt. 



Figure 3.1: Z is periodic over 
strips 



dz_ 

dt 



< 4a(0')0" - fa < 0. 



Therefore, at such a maximum point Z is non-increasing in time and so Z < 0. 

The reason for the restriction x ^ y is that is not differentiable here. When the 
maximum is attained at such a point, then Z(x, x, t) = 0, so in either case Z < for all 
t. The result follows. □ 

We can find explicit estimates for more general equations by choosing an explicit 
barrier. 



3.3 Description of a barrier 

This barrier, tp, will be used often. 

Let $ be the fundamental solution to the heat equation, 

so that & yy = 4c$ 4 , where c is a positive constant and t > 0. Implicitly define ip(z, t) by 

z = - l,t) - $0 + !,*)■ 
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This function has the property that as t -»■ 0, 

Y>0M) 



1 ^ > o, 

1 z < 0. 



and that ^(0, t) = for all t > 0. 
We can calculate 

^ = I 

$ y (V>- i,t)-^ v (V + i,t)' 

< = - o') 3 [* w (^ - 1, t) - <j> ra (v + 1, *)] , 

and (third derivatives are included for completeness but not used until a later chapter) 



//\2 



while 



' - 1, t) - + 1, *) " 



Routine calculations yield 



2c 



~ t 3 /2 

Ac 2 y 



2cy 
2c?/ 2 



cxp 



yyy 



t 5 / 2 
1 



-1 + 



2cy 2 



2cy 2 



t 

exp 
exp 
exp 



t 



(3.5) 



The partial differential equation satisfied by ijj is 

dip _ 1 i)" 
~dt ~ Ac^j/ 2 ' 



3.4 An explicit estimate for periodic equations 

Theorem 3.2. Let u : R x [0, T) ->■ R £>e a i? 2 solution of 

u t = a(u x ,u, x, t)u xx + 
n(-,0) = u 
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where u is continuous; both u and a are periodic, u (x + L) = u (x), a(p, q,x,t) = 
a(p, q,x + L, t) (and therefore u is also periodic); osc u(-,t) < M; and where we can 
find positive constants A and P such that 

a(p, q, x, t)p 2 > A for all \p\ > P. (3.6) 

Then there is a T > such that for t e (0, T'}, 

\u x \ < CiVt(l + t)exp(C 2 /t), 
where T', C\ and C 2 are dependent on M, A and P. 



Proof: Let Z be as before, with <j>(z,t) := 2Mip(z/2M,t/AM 2 ) for the V defined in 
Section 13.31 with the constant c given by c = max , CP 2 ) , where C will be chosen 
later. 

Consider the region 

G := { (x,y,t) : < y-x < z M (t),0 < t < 2cM 2 /3 } , 
where z M (t) satisfies 4>(z M (t),t) = M. Explicitly, 



ZM(t) = 




(3.7) 




Figure 3.2: The (periodic) region G at some time t 



As before, Z is periodic over strips parallel to y+x = and so it attains its maximum 
on G. We first show that Z < on the boundary of G. 

For y - x / 0, as t ->■ 0, <f> ->■ 2M and so Z < 0. 

At y - x = 0, <j)(0,t) = for all t > 0, and so Z(x,x,t) =0. At y - x = z M (t), 
(j)(z M ,t) = M and so Z < 0. 

Now suppose that Z attains a maximum on the interior of G. 
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It follows from (13.31 that at the maximum, 



dZ 

— < 2 [a (<f)',u(x,t),x,t) +a((j)',u(y,t),y,t)] (fi" - (j) t 



2 [a (ip', u(x, t),x,t) +a (ip', u(y, t),y, t)} 
1]/' 



V ipt 



2M 

The second derivative 



2M 2M 

2a (?/>', u(x, t),x, t) + 2a (t//, u(y, t), y, t) 



Acip l2 \ 



i>" = - (*l>'y [4>yy{iP - l,t/4M 2 ) - 4> yy (^ + M/4M 2 )] 

is negative, as ip' is positive, and the part inside the square brackets [ • ] is positive in 

G\U< 2cM 2 /3 = V. 



We can estimate 



ib' 2 > inf ib' 2 
G 



> 



( 

sup %(*p- l,t/4M 2 ) -$ y (^ + l,i/4M 2 ) 

i 0<^<l/2 
\ 0<t<T' 



>( sup sup |$ w (^- l,i/4M 2 )| + |$j,(V> + l,£/4M 2 ) 

\0<^<l/2 0<t<T' 



> sup 

\0<V><l/2 



Vo<vSl/2 (2/3)3/2^|^ - 1|2 ' (2/3)3/2^ + 1| 2 

2e 3 
P5 2 



*„U-l,-c(^-l) 5 



+ $ y U + l,-c(-0 + l) : 



2e -3/2 



+ 



2e -3/2 



> C 



where the last line follows by choosing C = 3 3 5 2 /(2e 3 ) and recalling that c > CP 2 . 



Now we can use the condition Q3.6| . controlling the degeneracy of a, to estimate 

1 



dZ ib" 
~dt ~ 2M 

~ 2M 
< 



2a (?// , u(x, t),x, t) + 2a (ip',u(y, t),y, t) 



4c0' 



/2 



ip' 2 4ap 



j 2 



as c > (16A)" 1 . So, Z t < at an internal maximum, Z < on the boundary, and so 
Z < on G. 



Chapter 3: Gradient estimates in one space dimension 



21 



Explicitly, for (x,y,t) e G, Z < means 

\y — x\ t 



u(x,t) - u(y,t) < 2M-0 



2M '4M 2 



= \y - x\ip' (0, 



t 

UP 

t 3/2 / 4c M 2 \ 

= |y " x| 4<2M)3 eXP l^J' ( ) 

which is an estimate for the difference quotient \u(x,t) - u(y,t)\/\x - y\. 

We can obtain an estimate for (x, y, t) G, y > x by observing that for z > z M (t), 

u(y, t) - u(x, t) - < u(y, t) - u(x, t) - M < 

ZM(t) 

so that as \y - x\ > zm(*)> 

u{y,t)-u{x,t) < < ly--l^exp(^). (3.9) 

So far, we have estimates for when y > x. We can find identical estimates for the 
region where y < x by reflecting in the line x — y = 0. 

Letting y -> x in the estimates for the difference quotients d3.8fr and i|3.9| gives the 
result. □ 

Comparing with the special barrier V> gives us a gradient estimate for solutions 
of quasilinear equations with this scaling. We will use this estimate in Chapter^] 

Corollary 3.3 (Gradient estimates for the barrier <p). Letip be a smooth solution of 

ft < a((ff,ip,z,t)(ff', 

on (0,oo) x (0,T), with the initial and boundary conditions 

<p(z,0) = 1 fort > andz / 0, 

<p(0,t) = fort > 0, 

ip(x,t) — > l as x — >• 00 fort > 0. 

/f 

a(p, r, sc, t)p 2 > A > for a// |p| > P > 0, 
then there is a T > sue/) fr?af for t e (0, T'] 

(f/{0, t) < Ci Vt(l + 1) exp (C 2 /t) 

i/i/r»ere X", C\ and C 2 depend on A and P. 

Proof: We apply the comparison principle to ip and ip € (x,t) ■— e(l+t) + 2ip(x/2,t/4), 
where ip is as in Section l3~3l with the constant c = (4A) _1 . 
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The barrier ip e dominates <p on (0, oo) x {0} and {0} x [0, T}. 

Choose T > and z\ so that on (o,*i) x (0,1*), ^" < 0, V e/ > P and at z x , 

4> e {zi,t) > l > 99(21, t). 

Then = -y? t + a(<//, 99, z, t)<p" > and 

Pip e = -e 



<0, 

so Theorem EUJ implies that ip e > 9? on (0,zi) x (0,T') for all e > 0, and as e ->• 0, 
■0° > <f. Since 0°(O, t) = (f(0, t) this gives us the boundary gradient estimate 

<p'(0,t) <^'(0,t/4) < CiVt(l + t)exp(C 2 /t). 

□ 



WO 



M2 



4c 



a(il}' ,il) e ,x,t) 



3.5 Interior estimates for non-periodic equations with 6 = 

Theorem 3.4. Let u: fix [0, T] -> R be a F 2 solution to 

u t = a(u x ,u,x,t)u xx , 

where OcRis an open interval and where there are positive constants A and P so 
that 

a(p,q,x,t)p 2 > A forall\p\ > P. (3.10) 
/f oscq x [ 0)T ] u < M, then we can find an estimate for < t < T 



u'(x,t)\ < CiVi(l + i)exp 




where T', C\ and C 2 are dependent on A, P, M and dist(x, dfl). 

We modify the previous proof, introducing a new boundary in the x,y coordinates, 
since we no longer have compactness of the domain through periodicity. We will seek 
to avoid a maximum of Z occurring on the new boundary. 

Proof: Firstly, suppose that = (-1, 1). 
We consider the sub-region 

G:={(x,y,t) € (-1,1) 2 x [0, T'\ : < t < T' < T, \y + x\ < 1, < y - x < z M (t) }, 

where z M (t) is as before in i|3.7| and V is chosen so that z M {t) < 1 for t < T . 
Define Z on G by 

Z(x, y, t) := u(y, t) - u(x, t) - <p(x, y, t). 

In order to avoid a positive maximum of Z on the boundary, we will ensure that 4> 
satisfies 
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• 4>(x, y, 0) > M for y ^ x 

• (f)(x,y,t) > M for \y — x\ = z M (t) 

• <j)(x,x,t) > for t > 

• <j)(x, y,t)> M for \y + x\ = 1. 



/G 


>\ / 

A 

\ 


X / 

\\ / 

V \' 













Figure 3.3: The region G at some time t 



We choose 



= 2MV> 



+ 7 (x + y-2P) 2 , 



2M ' 4M 2 

where ^ is the explicit barrier defined in Section 13.31 for positive constants c and 7. 
We will also choose with |/3| < 1/2 later. 

This 4> satisfies the first three conditions. In order to fulfill the final condition, choose 

7 = M/(l - 2\p\) 2 . Then at \x + y\ = 1, 



7 (x + y-2/3) 2 > (j(\x + y\ 



7(1 



M. 



Now suppose that Z first reaches a positive maximum at an internal point (x, y, t) e 
G. As usual, we calculate that at this point first derivatives are zero 

= Z x = -u'(x, t)+ i>' - 2 7 (x + y- 2(3), 
= Z y = u'(y, t)-i//- 2 7 (x + y - 2/3), 

and the matrix of second derivatives is negative semi-definite 

V(s,t)- 5^" -2 7 217^' 
^"-27 «"(y,t)-j 

We use these in the evolution equation for Z 

dZ 



> 



Zxx 


Zxy 




Zyx 


Zyy_ 





27 



at 



a(uy,u(y),y,t)u yy - a(ux,u(x),x,t)u xx - 4>t 



24 



a(uy,u(y),y,t) 



2M 



Z yy + 777 + 2 7 



a (iia.,u(:c),2;,t) 



— - (a(u y ,u(y),y,t) + a(u x , u(x) , x , t)) Z, 



■r.V 



+ (a(u y , u(y), j/, i) + a(« x , x, t)) I — - 27 



trace 



a(u x , u(x) , x , t) —a(u x ,u(x),x,i) 
-a{u y ,u(y),y,t) a{u y , u(y) , y , t) 



Z X X Z X y 

Zy X Z yy 



1p" 1p" Ipt 

+ 2a(u v , u{y). y. t)—— + 2a{u x ,u(x),x, t)—— — 



< 



2a(u y , u(y),y, t) + 2a(u x ,u(x),x, t) 



1 



4ctp" 



2M' 



the last line applying at the internal maximum . As before, 4>" < for t < T < 
2cM 2 /3. The first derivative condition for Z implies that u x = ij/ - 2^{x + y - 2(3) and 
u y = ip' + 2-y(x + y - 2(3) and SO 

max(\u x \,\u y \) >\ip'\> P, (3.11) 

where the final inequality comes from choosing c > CP 2 as in the previous section. 
We can then exploit i|3.10| . the condition on a; 

2A 



2a (u y , u(y),y, t) + 2a (u x ,u(x),x, t) > 



> 



max (\u x \,\u y \) 2 
2A 

(\i;'\+2 1 \x + y-2p\) 2 - 



This is greater than 1/Acifj' 2 whenever (\/8Ac- > 47, so we use (|3.1 T) , the lower 
bound on and choose 



c > 



(4 7 + Pf 
SAP 2 



AM 



(1- w 



+ p 



1 



SAP 2 ' 



Now Z t < at interior maxima, and so the parabolic maximum principle ensures 
that Z < on G. 

For points outside G, but in the rectangle \y + x\ < 1, z M {t) < y - x < 1, 

M|y — x\ 



u(y,t) - u{x,t) <M< 



z M (t) 



y/t fcM 2 ' 

< exp \ \y — x\ 

~ 2M \ t 

We can repeat both these estimates for a reflected region where x > y, putting them 
all together gives 



\u(x,t) - u(y,t)\ < 2Mtp 



, M(x-y) 2 y/t (cM 2 \. 

+ ^TTTTTTT + 7777 eX P ]\y - X \ 



( \y — x\ t 

V 2M ' 4M 2 J ' (1-2|/3|) 2 ' 2M 



t 
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Now, for \y\ < 1/2, set j3 = y and let x -> y to give a gradient estimate at y 

t 3 ' 2 (AcM 2 \ /cM 2 \ 

^ 24^ eX H^ J + 2M eXP ^ ) ■ 



For the case of a general interval O = [sci, rr 2 ], we can rescale around a point j/efl 
by using scaled coordinates x = 2{x — y)/dist(y, dSl). We obtain the estimate 



\u'(y,t)\ < 



dist(y, d£l) 

where c depends on ^4, P, M, and dist(y, eft}). □ 



t 3 ' 2 /AcM 2 \ „ r (cM 2 \ 

2w ex H^J +M ^ exp l^J 



3.6 A generalisation to fully nonlinear equations 



In this section we consider equations of the form 

ut = F(u xx ,u x ,u,x,t), (3.12) 
where F : R 3 x n x [0, T] -> R is C 1 . 

Let n be a smooth solution to Q3.12| . As before, define 

Z(x,y,t) := u(y,t) - u(x,t) - <j>(\y - x\,t), 



and suppose it first becomes non-negative at some point (x, y, t), with y > x. First and 
second derivatives of Z will satisfy i|3.1| and <3.2| . but the evolution equation for Z will 
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be given by 

^ = u t {y,t) - u t (x,t) - <j) t (\y - x\,t) 

= F(u yy ,Uy,u(y),y,t) - F(u xx ,u x ,u(x),x,t) - <p t (\y - x\,t) 
1 dF 

-j-(su yy + (1 - s)u xx , suy + (1 - s)u x , su(y) + (1 - s)u(x) , sy + (1 - s)x, t)ds 

- <M|y - z|,t) 

c] / + (1 - S)U XX , . . . ) ds 

f 1 dF 

+ [Uy ~U X ] J — (sUyy + (1 ~ s)^, . . . ) (fa 

+ - u(x)J / — (sUyy + (1 - s)^, • • • ) ds 
Jo a Q 

f 1 dF 

+ [y-x) J — (suyy + (1 - s)u xx , . . .) ds - <j)t(\y - x\,t) 

\ pi Qp 

j+2(f)"J — (suyy + (1 - s)u xx ,. . .) ds 
f 1 dF 

~ (Cl + C 2 )(p" - 4> t +4> -TT- (SUyy + (1 - S)U XX , ...)ds 

Jo a Q 

f 1 OF 

+ [y-x] J — {suyy + (1 - s)u xx , ...)ds, 

where we have used that u y = u x at a spatial maximum, and have abbreviated 

f 1 dF 

J — (suyy + (1 - s)u xx , . . .) ds = a 
and where we have added and subtracted (ci + c 2 )Z xy for some c x , c 2 . If we choose 



trace 



( 


a 


Cl 




Z xx 


Z x y 




C 2 


a 




_Zy x 


Zyy_ 



c\ = C2 = -a 



the first matrix above is positive semi-definite. Since the matrix of second derivatives 
is negative semi-definite, we have 



d 4<-^" rl9F 



(SUyy + (1 - S)U XX , ...)dS-(j>t 

o ar 



f 1 dF 

+ c t > ~w ( su yy + s ) u xx, ■ ■ • ) ds 
Jo °1 



f 1 dF 

+ ^ ~ X V !)x~ ( SUyy + ^ 



s)«xx, . . . ) ds. 



If we make quite harsh restrictions on F, then we can use our explicit barrier to find 
an analogue of Theorem I3T21 for periodic nonlinear equations. 



Theorem 3.5 (Nonlinear version of Theorem ET2l. Letu :Rx (0,T) ^ M. be aC 2 
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solution of 

u t = F(u xx ,u x ,u,t) 

where u is continuous and periodic, u (x + L) = u (x), (and therefore u is also peri- 
odic); osc u(-, t) < M; where we can find positive constants A and P such that 

dF 

-Tr{r, P ,q,t)p 2 > A forall\ P \ > P; 

or 

dF 

and where — < . 

oq 

Then there is a T > such that for t e (0, T'], 

\ux\ < CiVt(l + t)exp(C 2 /t), 
where T', C\ and C 2 are dependent on M, A and P. 



Proof: The proof is the same as that of Theorem l3T2l including the choice <j>(z,t) = 
2Mip{z/2M,t/AM 2 ), but instead of using inequality (|3.3| for interior maximum points, 
we have 



<o, 



.4 



y 2 4c</> /2 



where the omitted argument of the derivatives of F, denoted by (. . . ), is (su yy + (1 
s)u xx , suy + (1 - s)u x , su(y) + (1 - s)u(x),sy + (1 - s)x, t). □ 



Chapter 4 



An existence result for a parabolic 
equation in one space dimension 

Although this is a standard result (see Theorem 12.25 of (25|), for completeness we 
sketch a short time existence result in the one-dimensional case, where the spatial 
domain is £1 = (x , x%) c R, and the initial and boundary data is continuous. 
The parabolic equation is 

u t = a(u x ,u,x,t)u xx , (4.1) 
with initial and boundary data prescribed by 

u(x, t) = u (x, t) for (x, t) G P{tt x [0, T]). (4.2) 

We require that a > is in H a (K,) for all bounded /CcixRxfix [0,T] and some 
a G (0,1). 

This implies that for every such K, we can find positive A/e and such that 

\k < a(p, q, x, t) < Ajc, when (p, q, x, t) g JC. (4.3) 

When we can find bounds of this type that depend only on the gradient, we will 
write 

\{K) < a(p, q, x, t) < A(K), for \p\ < K. (4.4) 

Suppose also that there are positive constants A and P such that 

a(p, q, x, t)p 2 > A > 0, for \p\ > P. (4.5) 

The first part of this chapter is a survey of the main steps needed to find the exis- 
tence result for the Cauchy-Dirichlet problem with H 1+/3 initial and boundary data. We 
follow the treatment in Lieberman |25] . 

These results mean that when we approximate continuous initial data by smooth 
initial data, a solution will exist for the approximate initial data. In the later parts of the 
chapter, we use the gradient estimate established in Chapter|3]to find uniform gradient 
estimates for t > 0. This will gives us a solution for t > 0; in order to show that this 
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approaches the initial data as t -> 0, we will need some displacement estimates which 
limit the distance a function can travel in a given time. 

4.1 Existence of solutions with //, , initial and boundary 
data 

Theorem 4.1 . Consider the Cauchy-Dirichlet problem given by \4A) and \4.2) . where 

<p< 1. 

Suppose that u is defined on the parabolic boundary V (0 x [0, T\) and u e 
Hi + p(T>). Also, suppose that either u is time-independent, or else there are con- 
stants A and P such that \4.5) is satisfied. 

Then there is a smooth solution u e C 2+1 (fi x (0, T)) n C (H x [0, T]) . 

Tn/'s solution has a gradient bound \u\\ + s,s/2 < C where C depends on \ u \ 1+f3 ^/ 2 , 
[3, Xfc, Ajc and diam Q. 

The proof of this result follows a standard pattern for showing existence — a bound 
on sup \u\ ; a bound on sup \Du\ ; a Holder gradient bound \Du\ a ; and then the applica- 
tion of a fixed point theorem. These steps are sketched by the following results. 

We begin by using the comparison principle to bound \u\. 

Lemma 4.2 (A bound on sup \u\). Ifu is a smooth solution of J4J), \4.2\ in O, x [0, T], 



Proof idea: Set k = supn + and apply the comparison principle (Theorem |2jJ to 

u and k on E := { (x, t) e!Jx [0, T] : u(x, t) > }. Since k > u on dE, it follows that 
k > u on all of E and so on all of Q x [0, T\. 

Similar steps can be followed to find that inf^ X [ 0T ]n > infu _ , completing the 
result. □ 

We begin our gradient estimates with a boundary gradient estimate. 

Lemma 4.3 (Boundary gradient estimate). LetO < (3 < 1. If u is a smooth solution 
of (|4.1| . \4.2) in O, x [0,T], and either u is time-independent or else a satisfies the 
condition \4.5) . then 



In fact, we can relax the regularity requirements on the initial and boundary data 
and still find a continuity estimate on the boundary. 

A modulus of continuity is a concave, continuous function u : R + — >■ R + , with 
cj(0) = 0. This uj is a modulus of continuity for a function g at y if 



then 



sup t)\ < sup |uo 

Qx[0,T] 




)l 



< L 



g{x) - g(y)\ < u(\x - y\) 



for all x in the domain of g. It is a modulus of continuity for g if the above relationship 
also holds for all y in the domain of g. 
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A modulus of continuity can be defined for every continuous function on a closed 
bounded set. 

Lemma 4.4 (Boundary continuity estimate). Let u satisfy <4.1| . (|4.2| . where u has 
modulus of continuity u, and suppose there are positive constants n and P so that 

\p\k{p, q,x,t) + 1 < na(p, q, x, t)p 2 (4.6) 

whenever \p\ > P. 

Then u has a modulus of continuity on the boundary 

\u(x,t) - u(y,t)\ < uj*(\x - y\) 

forx e andy e dQ, where lo* can be determined byco, sup-p \u \,fl, and a. 

Equipped with the boundary gradient estimate, we can now find a global gradient 
estimate. In this one-dimensional case, the global gradient estimate is the result of 
Kruzkov mentioned in Chapter|3l 

Lemma 4.5 (Global gradient estimate). If u is a smooth solution of {HP , (|4.2| in 

Qx[o,T] with an oscillation bound osc u = M, and a Lipschitz estimate on the parabolic 
boundary \u(x,t) - u(y,t)\ < L\x - y\ for all (x,t) e Vip, x [0,T]) andy e £1, then 

sup |ti x | < 2L. 
nx[o,T] 

Lemma 4.6 (Global Holder gradient estimate). Suppose thatu satisfies Q4.1| . \4.2\ 

in U x [0,T], where there are positive constants and A/c such that whenever 

(p, q, x, t) in the set K •— { (p, q, x, t) : \p\ < K, \q\ < M, x e O, t e [0, T] }, 

A/c < a(p,q,x,t) < A K . 

Ifu g C 2+l {n x [0,T]) n C(Tl x [0,T]), se? M = sup|n| andK = sup|D«|. Then 
there are positive constants a and C determined by (3, \ K , K K and diam Q, such that 

[Du] a < C (sup \u\ + sup \Du\ + |lio|l+/3,/3/2) • 

Now that we have bounds for \u\ 1+a ^ a / 2 , we can apply the following existence the- 
orem, which is derived from a fixed point theorem. 

Lemma 4.7 (Existence theorem). Letu be in H± + s for some 5 e (0, 1) 

If there is a constant M s independent ofe such that any solution of Q4.1| . Q4.2| on 
Qx [0, e) satisfies 

l«|l+<5,5/2 < M s, 

then there is a solution of the Cauchy-Dirichlet problem (|4jJ, \4.2) in U x [0, T}. 

4.2 Displacement estimates 

The following estimates for the displacement suffered in a given interval of time by a 
function moving under a parabolic flow apply to any strictly parabolic operator satisfy- 
ing bounds of the form (|4.3| or i|4.4} . 
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Lemma 4.8 (Displacement estimate for Lipschitz initial data). Letu -. £lx[0,T] — ► R 
satisfy \4A) . where fici, and a has bounds of the form (|4.4f . 

Suppose that u has initial data whose graph lies below a cone centred at some 
point h 

u(x, 0) < L\x — h\ 

and, in the case that tl^R, whose boundary data lies below the same cone 

u(x,t) < L\x — h\, x £ dtt. 

Then, at later times, 

u(x, t) < L(x - h)£ri ( + 2L 

where A = A(L) is given by (|4.4| . 
Proof: For some small e > 0, set 

v(x, t) := L(x - h)£rf (J^(x - h)\ + ^/^xp L&^tj , 

which satisfies the heat equation 

Vt = — C 

and approaches the cone of gradient L centred at h as t + e — >• 0. 



Figure 4.1 : The barrier v(-,t) 



Note that u(x,0) > L\x - h\, that \v x (x,t)\ = L Sri (xyjc/{t + e)) 
v xx > 0, SO 



< L and that 



v t - a(v x , v, x, t)v xx >v t - sup a(p, v, x, t)v x 

\p\<l 

>v t - A(L)v xx 

= 
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where we choose c _1 = 4A(L). 

The estimate follows by applying the comparison principle (Theorem l2~2fr to show 
that u(x, t) < v(x, t), and then letting e ->■ 0. 

□ 

Now, we apply this to three different cases, firstly when u initially satisfies a Holder 
condition and when we have polynomial growth in A, secondly when u(-,0) has a 
modulus of continuity, and thirdly when u is initially bounded by a step function. 



Corollary 4.9 (Displacement estimate for Holder initial data). Let u:Kx[0,T] — > R 
satisfy S3). Suppose that a not only satisfies \4A\ . but more specifically has at most 
polynomial growth, so that 

a{p,q,x,t)<A(l + K m )for\p\<K (4.8) 

where A and m are positive constants. Also, suppose that u(-,0) satisfies a Holder 
condition around some point h 

KM) -u(x,0)\ < L\x- h\ a , < a < 1. 

Then, af /ater times, 

\u(h,0) - u(h,t)\ < c(a,m,L, A)t 2+" 1 " 1 -") . 



Proof: For simplicity, assume = and u(h,0) = 0. The initial data is bounded 




Figure 4.2: The bounding cusp is itself bounded above by cones 



above by cones centred at h and indexed by k, the (positive) x-coordinate of the point 
of contact with the bounding cusp L\x\ a , so 

u{x, 0) < L\x\ a < aLk a ~ l \x\ + L(l - a)k a . 
The estimate taken at x = is then 

:a-l I -i , /r7.a-l^ m ^ 1 /2 , / At 



«(0,t) < 2Lak a ~ (1 + {L^-'y ) ' \j — + L(l - a)k a 
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and optimizing over k gives 

u(0,t) < c(q, m, L, K)t 2 + m " 1 - a ) . 

□ 

Corollary 4.10 (Displacement estimate for continuous initial data). Suppose that 
u : R x [0, T] — > E satisfies (Q) and <4.4| . wnere u nas /M/'a/ dafa i/w'fn a modulus of 
continuity to at a point h 

\u(h,0) - u(x,0)\ <uj(\x- h\) . 

Then 

\u(h,0) - u(h,t)\ < c(t) 
where c is dependent on oj and A, and where c(t) — > as t — >• 0. 

Proof: For simplicity, assume h = and n(/i, 0) = 0. Consider the cones 

Cfc(x) = c fc (\x\ - k) + 

indexed by > 0, the (positive) x-coordinate of a point of contact with u. As u is 
concave it has both left and right derivatives, and we can choose the slope of the cone 
cfe = J_{k). Then 

u(x,0) = u(x,0) - u(0,0) < uj(\x\) < u'_(k) {\x\ - k) + u{k) = C k (x). 

Now we have a cone as an upper boundary, we can use estimate (|4.7| at x = 

u(0,i) < 2uj'_(k)]J^- - u'_(k)k + u(k), 

where A = A(wi.(A;)) is given by (|4.4| . Minimize this over A: to get the displacement 
bound 



c(t) := inf \2u/(k)\ —-uL 
k>o\ V vr 



In order to show that c(rj) as t -> 0, let 5 > 0. Asw is concave and positive, it 
has positive left derivative and for k > we have 

< oj'_{k)k < uj{k). 

And as uj is continuous, 

< lim u'_(k)k < lim u(k) = 0. 
Choose k = k s so that u(k s ) - uj'_{k & )k & < 5. Choose r so that 



6y/A{u/_(k S )) 
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then for all t <t, 



c{t) < 2u'_(k s )J^ - uj'_{k 5 )k s + (j(k s ) < 25, 



and so c(t) -> 0. □ 

Set cr to be the maximal monotone graph 



cr(x) = < 



+ 1, x > 

[-1,1], x = (4.9) 

-1, X<0 



which we will refer to as the step "function". 

Corollary 4.11 (Displacement estimate for step functions). If u satisfies <4.1| and 
\4A\ . and is initially bounded by a step function 



u(x, 0) < ca(x), 

then forx < 

, . ( 4e [M } 
u(x, t) < mm < - — - y c, c > , 



w/7ere A = A(2c/|x|) as /'n i|4.4| . 

Proof: Near some point h < 0, 0) satisfies a Lipschitz condition 

u(x,0) < L h |x-/i| -c 

where L fc = 2c/\h\. 




Figure 4.3: Cone bounding the step function 



Chapter 4: Existence for one space dimension 



35 



LemmaEiHIthen gives that 



(x, t) < inf — 
h<0 \h\ 




x — h 



) 



4c [At 



(x - hf 



) 



RVV exp 



— c 



and if we let h = x then we find that for x < 0, 




— c. 




□ 



4.3 Existence of solutions with continuous initial data 

Theorem 4.12. Consider the Cauchy-Dirichlet problem given by \AA) and \A.2) . If 

u e C (V(£l x [0, T])) and if there are constants A and P such that (|4.5| holds, then 
H?n . S3 ^as a solution u e C 2+1 (n x (0, T)) n C (H x [0, T]) . 

The first step in the proof of the above is to approximate u by u € in C°°, so that 
su Pxen \ u o~ u o\ < e. 

Lemma 4.13 (Existence of solutions with approximate boundary data). For all 

e > 0, there exist solutions u e -. Q x [0, T] — > R to jQ) vv/to boundary data uq These 
solutions are in C 2+1 (Q x [o, T]) n C (H x [0, T]) . 

Proof: As u e Q is in the Holder space H 1+/3 , this is a consequence of TheoremEUJ □ 
Lemma 4.14 (Existence of uniform oscillation bound). For all e > 0, 



Proof: For any fixed e, set k = sup u e + and apply the comparison principle (Theorem 

HUJ to k and u e on E = { (x,t) e O x [0,T] : u e {x,t) > 0}. Since jfe > u e on it 
follows that k > u e on all of E 1 and hence on all of CI x [0, T], and so 



where the last inequality will hold for small enough e. 

This leads to a uniform oscillation bound for u e , which we denote by M — 



oscm 6 < 4 (sup | i*o I ) - 



sup \u e (x,t)\ <sup|uq| <2(sup|uo|) 

f2x[0,T] 



oscu e < 4 sup I no I =: M. 



□ 



36 



Theorem |331 gives a uniform gradient bound on interior sets, up to some time 

T > 0. For t e (0,772), 



where T', C\ and C 2 are dependent on A, P, M and dist(0',90). 

Lemma 4.15 (Higher regularity on interior sets). On interior sets W x (2t ,T ! ) we 

can estimate higher derivatives 

\u t \2+k+a < C 

where C depends on dist(fi', dQ), diam(fi), t , A, P, \a\ a andM. 

Proof: Once we have an oscillation bound M and a gradient bound L(t Q ), <4.3| im- 
plies uniform parabolicity. A uniform Holder gradient bound on interior sets results 
from Theorem 12.2 of |25]. In particular, on interior sets and when T'/2 > t > 0, 



and C also depends on L(t ) + M, and diamQ. 

Equipped with a Holder gradient bound, we can treat the equation as a uniformly 
parabolic equation with Holder continuous coefficients, and use standard results, such 
as Theorem lA.4l to find that u f - is uniformly bounded in H 2+a x {2to,T')). 

From here, it is possible to use the bootstrapping method to obtain interior esti- 
mates for all higher derivatives. 



Corollary 4.16. On any interior set n'x(f , T'), there exists a subsequence converging 
to some u that also solves the partial differential equation Q4.1| . 

In order for this u to be a solution of the Cauchy-Dirichlet problem, we need to 
show that u attains the initial and boundary data. 

Lemma 4.17 (Convergence to initial data). On any spatially interior set fi\ 




[u e x]a-,n'x{2to,T') < Cmm{dist(Q',dtt),y/t^y 
where both a and C depend on \ K and A^, given by (43J, with 



/C = { (p, q, x, t) : \p\ < L(t ), \q\ < M, x E Q, and t E [0, T] } 



□ 



sup \u(x, t) — Uq(x)\ — > as t — > 0. 
xen' 



Proof: Let x be any point in O'. Let u be a modulus of continuity for u . 
We can off-set u e by defining 



w e {y,t) := u e (y,t) 



u e (x,0) + «o(x) 
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so that w e (x,0) = u (x). Let u be the limit of a subsequence u e , as in Corollary 14.161 

\u(x, t) — uq{x)\ = lim \u e (x, t) — uq{x)\ 

= lim \ w e (x, t) + u e (x, 0) — 2uq(x)\ 

< lim (\w e (x, t) — uq(x)\ + \u e (x, 0) — uq{x)\) 

= lim \w e (x, t) — w e (x, 0)| + lim \u 6 (x, 0) — uo(x)\. 

e— >0 e^O 

The second of these terms is zero. To estimate the first term, note that the approxi- 
mations u e (-,0) satisfy the same the same continuity condition as n , and therefore so 
does w e (-, 0), with |u> e (0, x) - w e (0, y)\ < uj{\x- y\), for all x, y £ Q. Corollarv R.IOI then 
gives the estimate 

\w e {x,t) - w e (x,0)\ < c(t), 



where c depends only on the exact forms of u and A/c (given in (|4.3|). In particular, c 
is independent of x and e, and as c(t) -> as t ->■ 0, the result follows. □ 

More specific continuity-in-time estimates are given by the continuity of the initial 
data and the upper growth bound of a. If, for example, the initial data is Holder contin- 
uous 

\u (x) - u {y)\ < L\x - y\ a 



and a has polynomial growth in the gradient term, satisfying 0~8) for constants A and 
m, then Corollary l4~9l indicates that c(t) = C(a, m, L, A)\t\ 2+™a-°o . 



Lemma 4.1 8 (Convergence to boundary data). We can continuously extend u(-,t), 
defined on the interior of£l at time t, to Tl. Moreover, u = u on the boundary. 

Proof: We need to show that for y e dO,, lim^j, u(x, t) = u (y, t). 
We note that our parabolic equation satisfies condition since 



2 2 

MA(p, q, x) + 1 = \p\\a(p,q,x)\ + 1 < —a\p\ 



for \p\ > P, using (475). 

Let u> be a modulus of continuity for u . As each u e has at least the same modulus 
of continuity as n , Lemma l4~4l gives us an estimate uniform in t and e, 

\u e (x,t) -u e (y,t)\ <u*(\x-y\). 

Then for a point y e d£l and fixed t, 

sup \u(x,t) - u (y,t)\ = sup | \imu e (x,t) -u e (y,t) +u e (y,t) -u (y,t)\ 
B r {y)nn B r (y)r\Q 

< sup u*(\x-y\) 
B r (y)nn 

= w*(r) 

so as |x — y| — >• 0, oj*(\x-y\) — > and u(x,t) n (y, — that is, we can continuously 
extend u to u on dtt for t > 0. □ 
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4.4 Existence of entire solutions with stepped initial condi- 
tions 



Consider equation j4T), under the conditions on a given by Q4.3| and (|4.5| . 

Lemma 4.1 9. There exist entire solutions to this equation with the periodic, crenellated 
initial data 

g R (x) = Ma (sin(7rx/ J R)) , 



where a is given by Q4.9| . 



Proof: If we let g e be the smooth mollification of g R , then for |x| < R, 

Ma(x - e) < g e (x) < Ma(x + e). 

Theorem l4.1l ensures that there is a smooth solution u e to !4.1 I with initial condition 
g e , with a Holder gradient bound dependent on e. The gradient bound in Theorem l3~2l 
is independent of e; for t e (0, T), 

\u%\ < C\\ft (1 + t) exp(C 2 /t) 

where T , d and C 2 are dependent on M, A and P, but not fi. Higher gradient 
bounds for t > follow from the interior estimate i|4.15^ and we can find a subsequence 
converging to il r which also solves the equation on [t, T']. 

To show convergence of u R to the initial data, suppose that -R/2 < x < 0. As in 
Section l4~2l we can bound the initial data g e by cones centred at h e {-R/2, -e) — 

M 

9 e (x) < t, ;\x -h\- M. 

K ' ~ \h + e| ' ' 

Applying LemmaEiUto this, and setting h = x, we find that for -R/2 < x < -e, 

\x + e V vr 



where A = A (M/|x + e|) as in ^4.3| and so we have the estimate 

\u R (x, t) - g R (x)\ = | lim u e (x, t) + M\ 

2M /At 
|x + el V vr 
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A similar estimate holds for all x / nR, and so for all fi > 0, we can find t (dependent 
on x) such that \u R (x,t) - gn(x)\ < \i. □ 

Corollary 4.20. There exists an entire solution to this problem with the initial data 

uo(x) = Ma(x). 
This solution has a gradient estimate fort < T': 

K| < CiVt{l + t)exp(C 2 /t), 
where T', C\ and C 2 are dependent on M, A and P. 

Proof: Take the limit of the solutions u R given by the previous lemma as R — ► oo. □ 



Chapter 5 



Gradient estimates for parabolic 
equations in higher dimensions 



In this chapter we extend the methods of Chapter|3]to higher dimensions. 
Consider a smooth solution u : R n x [0, T] — > R to 

u t = a ij (Du, t)DijU + b(Du, t), (5.1 ) 

where A(p,t) = [a ij (p, t)] is a symmetric, positive semi-definite n x n matrix that is 
smoothly dependent on (p,t) g R n x [0,T]. 

Define 

a(p,t):=N 2 inf . (5.2) 



Compare this definition to that of the Bernstein £ function, (see Chapter 1 of [1 6]) 

£(p,q,x,t) = a tJ (p,q,x,t)piPj. 
Clearly, a(p)\p\ 2 < £ and if A, A are the smallest and largest eigenvalues of A, then for 

A < o(p) < ^2 < A- 

The middle inequality here becomes an equality when p is an eigenvector of A. 

Our aim is to reduce the n-dimensional problem to a parabolic equation in one 
space dimension; we can do this if a(p) is bounded below by a positive function of \p\. 
We will call this 

a(s) := Jnf a(p)- (5.3) 

:\p\=s 



For the existence of specific barriers we will require a control on the degeneracy of 
A — the existence of positive constants A and P such that 

a(s)s 2 > A for s > P. (5.4) 
Example 1 : If p is an eigenvector of A(p), then a(p) is the associated eigenvalue. 
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Example 2: As a specific example of the above, if a ij is of the form 

a ij (p) = «oo(p) (V j -yf)+ oo(P)S ( 5 - 5 ) 

for functions aoo, a : R n — > R, with a > and > 0, then a(p) = a (p). 
In the mean curvature flow case, = 1 and 

a (p) = a(p) 



1+ P 



Example 3: In the most general situation, if are the non-null eigenvectors of 
A(p) with eigenvalues A^(p) > 0, then 



a(p) 



\p\ 2 (E k {jJ jf-) 1 ifpe(NulL4(p))- 
otherwise. 



Example 4: In the case that A is positive definite, all eigenvalues are positive and 

M N 2 



As Example 2 shows, A need not be positive definite. 

Example 5: An elliptic operator [a ij ] is called of mean curvature type if there are 
positive constants A, A so that 



where m lJ are the coefficients of mean curvature flow Pbll271fT6| . For such equations, 
one can (under some conditions, particularly on the shape of the boundary) find apriori 
estimates on \Du\ in terms of \u\. It is therefore interesting to note that if [a tj ] satisfies 
only the lower inequality above, then it also satisfies j5.4| . 

Example 6: If one may be forgiven for referring to a future section, note that if 
the flow is the anisotropic mean curvature flow 1)9.5) of Section IPl then Lemma l931 
implies that S(|p|) > A \p\ 2 when F(p) > P. The function F will be positive and 
homogeneous of order one, so that ci|p| < F(p) < c 2 \p\. 

Thus, anisotropic mean curvature flow satisfies conditions i|5.3| and \5A) . 



5.1 Reduction to a one-dimensional problem 

Let u : R n x (0, T] — >■ R satisfy (Q), where A = [a^] is a symmetric, positive semi- 
definite matrix with a > 0. 

In the following, we generalise the calculations of Section l3~Tl to higher dimensions. 
As in the one-dimensional case, we begin our discussion by defining 

Z(x,y,t) := u(y,t) - u(x,t) - <p(\y - x\,t), 

where <p : R x [0, T] — >■ R is a C 2 function that will be chosen later. 
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At an internal maximum point of Z, the first derivative conditions are 

,y l — x % 



D xl Z 



DyiZ 



-Diu(x, t) — D x i(p(\y — x\),t) = —Diu(x, t) + <f> 
Diu(y,t) - D yi (/)(\y - x\,t) = Diu(y,t) - ft 



\y ~ x\ 



y — x 
\y - x\ 



0. 



(5.6) 



The second derivatives of Z are 

D xi x i z = -Uij(x,t) - D x i xj (j)(\y - x\,t) 

= -u tJ (x,t)- ^ yl -^- Xl) 
D y i yj Z = Uij(y,t) - Dy iy j(j){\y - x I , t ) 



= uij(y,t)-cf>' \ y - x \2 

D x i y jZ = —D x % y j(j){\y — x\,t) 

Jtf - x l ){y'J - x'J) , <j>' 



\y - x 



' (y l - x l ){yi - X?) 

Oin 



,,{y l - x l ){y 3 - x 3 ) 0' 



\y-x\ 



\y — x\ 2 
(y l — x l )(yi — x 3 ) 



\y - x\ 



+ 



(y i - x l )(y j -x j ) 



\y — x\ 2 ' \y — x\ \ 3 \y — x\ 2 

and at a maximum point, the matrix [D 2 Z] must be negative semi-definite. 



dZ_ 
dt 



The evolution equation for Z is 

ut(y,t) - ut(x,t) - <j) t 

a 13 (DyU, t)uij (y,t)+ b(D y u, t) - a lj (D x u, t)uij (x, t) - b(D x u, t) - <j> t 
a ij {D y u,t) [D yiyj Z + D yiyj ct>(\y - x\,t)] + b{D y u,t) 

- a lj {D x u,t) [-D x i xj Z- D x i xj (j)(\y - x\,t)] - b(D x u,t) 
-<k + 2c ij D xiyj Z+ 2c ij D xiyj (j)(\y - x\,t) 

(y l — x l ){y 3 — x 3 ) 



a lJ (D y u,t) 



D y i y j Z + : r= h i bij 

\y — x \ \y — x \ 



\y - x\ 



a ij (D x u,t) 



„(y i -x l )(yi -x j ) 
6 i 12 

\y — x\ z 

(y { - x { )(y j - x 3 ) 



\y-x\ V %3 \y - A 2 

+ b(D y u, t) - b(D x u, t)-(j)t + 2c ij D xiyj Z 



+ 2c 13 



^ {y l - x l )(y 3 - x 3 ) _ 4>' _ (y l - x % ){y 3 - x 3 ) 

\y — x\ 2 \y — x\ \ 13 \y — x\ 2 



where we add and subtract cross derivative terms with yet-to-be-chosen coefficients 
c tj . If we write £ = (y—x)/\y—x\, and assume that we are at an internal maximum, then 
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i)5.6> implies that Diu(x, t) = = Diu(y, t), and we can continue the calculation: 

dZ 

— = a l V£, t)D yiyj Z + o«(^'e, /)/>, , Z + l&D^Z 



\y — x 

+ b{(p'i,t)-b{(p'i,t)-4> t 

'A(4>'t) c 



(a^)-^) 



trace 



+ 2, 



|y - a: | 



d 2 ^) + 20" (^ T ^ - e r ce) 

(traced - fA£ - traceC + £ T C£) - fa 



The idea now is to choose the off-diagonal block C = [c ij ] in such a way that the 
2n x 2n matrix 

'A C 
C T A 



A' 



is positive semi-definite, leaves the coefficient of 4>" positive and sets the coefficient of 

(j)'/\y - x\ to zero. 

The first and third of these requirements imply that C is given by c ij = a ij (4>'{,, t) - 
c£i£j for some o 0. We can check that 



£ T A£ - fC£ = c> 



and that 



trace C - £ T C£ = trace [a y ' - c - £j{a lj - c&£j)£j 



= traced - f A£. 

If we set c = 2a(<f)'t;,t), with a defined by \5.2) . this maximizes the coefficient of 
b" , while keeping A' positive semi-definite. For any v,w e R n , 



(v T ,w T )A' ( v 



= v T Av + w T Aw + 2v T Cw 

= v T Av + w T Aw + 2v T Aw - 4a(£ • • w) 

= (v + w) T A(v + w) — a [(£ • v + w;) 2 — (£ • u — w;) 2 ] 

> (t> + ui) T A(?; + to) — a(£ • t> + w;) 2 

> 0. 



At the maximum point of Z, we find that 

In this way, we have reduced our problem to finding <p that satisfies the above 
equation, or, if we can find a lower bound on a dependent only on \p\, as in {53}, then 



45(^V' - <k < 0, 
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We will use the results of Chapter|3]to do this. 



5.2 Estimates for periodic solutions 



The following theorem is an analogue of Theorem 13721 for higher dimensions. In the 
special case of mean curvature flow, this is joint work with Ben Andrews. 
Let a be defined by \5.2) . 



a ij (Du,t)uij + b(Du,i) 



Theorem 5.1 . Let u:R n x [0, T] — > R be a smooth solution to 

du 

dt 

u(x, 0) = uq(x) 



where u is smooth with oscillation bound osc u < M, and is also spatially periodic, 
uq(x) =u (x + T), for some lattice T. 
Suppose thata(p) = a(\p\) > for all p. 

Iftp-.R + x [0, T] is a smooth solution to the auxilliary one-dimensional equation 

<Pt = 4a(WW, (5-7) 
and satisfies the boundary conditions 

lim ip(z, t) = 1 for z ^ 

(p(0, t) = for allt > (5.8) 
lim <p(z,t) -> 1 for allt > 

2— >00 

then 

\u(y,t) - u(x,t)\ < Mip ' 



M ' M 2 

Corollary 5.2. If there are positive constants A and P so that 

a(\p\)\ P \ 2 > A for\ P \ >P, (5.9) 
then there is a T > such that for t e (0, T'}, 

\Du\ < dVi(l + t) exp(C 2 /t), 
where T', C\ and C 2 depend on n, M, A , and P. 

Proof of Theorem Hj] The proof of this is substantially the same as the proof of 
Theorem l3~2l the gradient estimate for periodic, one-dimensional equations. 
As in the previous pages, let 

Z(x,y,t) := u(y,t) - u(x,t) - <f>(\y - x\,t), 

and choose (j>(z,t) = Mip (z/M,t/M 2 ), so that Z(x,y,0) < 0. 
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As u is periodic over the lattice r = {L u . . . , L n ), Z is periodic over regions 

{ (x, y, t) E K 2n x [0, T] : 2nL t - x t < y t < 2(n + l)U - x { }. 

On any one of these regions, note that Z(y, y, t) = and that Z(x, y, t) < M - 4> — ► 
as |yj - Xi\ — >• oo, so Z attains a spatial maximum on the region (and hence on the 
entire domain M 2n ). 

If there is a maximum point (x, y) at some t g (0, T') with x / y, then at this point 
Z is smooth and 

^ < i&^Of -4>t< 4fi(|y/|)g - g = 0. 

If x = y at the maximum point, then here Z(x, x, t) = and in either case, Z < 0. The 
estimate for \u(y, t) - u(x, t)\ follows. □ 

Proof of Corollary I5.2t If 5 satisfies the degeneracy condition ^5.9| . then for small 
times the gradient of <p may be estimated by Corollary l3T3l Letting x — >• y gives that 

\Du(y,t)\ < V(0,t) 

< CiVt(l + t) exp(C 2 /t). 

□ 

5.3 Estimates for boundary value problems 

In the special case of mean curvature flow the following theorem is joint work with Ben 
Andrews. 

Theorem 5.3 (Neumann problem). Let QcTtea convex domain with C 2 bound- 
ary, and let u be a smooth solution of 

dii ■ ■ 

— = a %3 (Du)uij + b(Du,t) 
D v u(x,t) = 0, forx £ dn,t > 0. 

If a ij and ip satisfy the same conditions as in Theorem \5l\ then for any x and y in Tt, 

\u{y,t) - u(x,t)\ < <p(\y - x\,t), 

where osc u = M. 

Furthermore if a lj satisfies the degeneracy condition (|5.4| . then fort e (0,T') a 
short-time gradient bound holds: 

\Du{x,t)\ < dVi(l + t)exp(C2/t) for(x,t) eSlx (0,T'] (5.10) 

where T , C± andC 2 depend on n, M, A , and P. 

Proof: As in the previous proof, set 

Z := u(y, t) - u(x, t) - M<p (\y - x\/M, t/M 2 ) . 
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Note that Z < when t = 0. 

For any t > 0, suppose that (x,y) is a spatial maximum of Z. We will consider 
the possibility that x and y are both interior points, that y is a boundary point and so 
is x, or that y is a boundary point while x is not (the converse follows without loss of 
generality). 

If both x and y are interior points, then the arguments of Theorem |5J] apply and 
Z < at this point. 

Consider the case that y is on the boundary dfl. If we take derivatives at y that 
are in directions // B that have no component normal to the boundary, then as before 
D fiy Z(x,y,t) = 0. On the other hand, let u y be the outward unit normal at y. The 
outwards-pointing derivative of Z here is 



d . , 
— Z(x,y + siy y ,t) 



= v y ■ Du{y,t) - ip 

s =o \y- x \ 
n / y- x 

= — ip 



\y - x\ 
<o, 

where we have used the boundary condition D Uy u(y,t) = and that as O is convex, 

(y — x) ■ v > 0. 

This inequality cannot be strict, for if it is, then there is a small s > such that 

Z(x,y- su y ,t) > Z(x,y,t) 
which would contradict that (x,y) is a maximum of Z. Therefore 



—Z(x,y + su y ,t) 



= 0, 

s=0 



and indeed D y Z(x,y,t) = 0. 

Now consider the position of x. If it is on the boundary, let v x be the outward unit 
normal at x, and so 



d 

— Z(x + sv x ,y,t) 



. y — x 

= —v x ■ Du(x, t) + ip ■ v x < 0, 

s=0 \y- x \ 



Again, this inequality cannot be strict if (x, y) is to be a maximum of Z, so the outward 
derivative ~rZ(x + su x ,y, t) = 0. As before, other non-normal derivatives are also 

zero, so D x Z(x,y,t) = 0. 

So, when both x and y are boundary points, DZ = and [D 2 Z] is negative semi- 
definite. We can argue as before that Z t < 0. 

In the case that x is an interior point, D X Z = and so DZ = 0, [D 2 Z] is negative 
semi-definite here, and Z t < 0. 

It follows that Z < for all t > 0. □ 

The highly geometric nature of mean curvature flow allows us to relax the condi- 
tions on the convexity of the boundary. In the following theorem we consider domains 
that are merely mean-convex. This means that at every point on the boundary, the 
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sum of the principal curvatures of 8ft is positive: 

n-l 

> 0. 

i=l 

Under the assumption of convexity (rather than mean-convexity), the following the- 
orem is joint work with Ben Andrews. 

Theorem 5.4 (Dirichlet problem for mean curvature flow). Let cl" be a mean- 
convex domain with a C 2 boundary, and letu be a smooth solution of the mean curva- 
ture flow for graphs 

du ( DiuDjU \ 

with prescribed boundary values 

u(x,t)=0 for x G d£l, t > 0, 
u(-,0) = ito. 

If ip is a smooth solution to curve shortening flow 

with boundary conditions given by <5.8| . then there is an estimate 

\u(y,t) -u(x,t)\ < 2M<p (k^,JL_) , 

where M = sup \u \. 

Proof: We find a boundary gradient estimate by defining a new Z B on Q x (0,T) 
which incorporates the distance to the boundary 

Z B (y,t) : =„( y ,t)-M^(A , 

where d(y) = dist(y, dQ) is a C 2 function in the neighbourhood of the boundary 

Q\n R := {y £Q : dist(y, dty < R }. 

Here, R = (sup an Kj) _1 and K t are the principal curvatures of <91X This close to the 
boundary, each point y has a unique closest point x e dQ. 

Choose T > so that for < t < T, if < d < Rthen ip' (d/(2M),t/(4M 2 )) > 0, 
and if d> Rthen cp(d/(2M),t/(4M 2 )) > 1. 

At t = 0, Z B < for all yeO. For y on the boundary, Z B (y, t) = - 2Mip < 0. 
For t < T and points at least distance R from the boundary, y e £1 R , Z B (y,t) < 
u{y,t) -M < 0. 
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We can find spatial derivatives for Z B : 

DyiZs = DyiU — cp' Did 
D ylyJ Z B = D y i yj u - -^DidDjd - cp Dijd. 



Now suppose that y is an interior maximum of Z B at some time t < T . At this 
point, DZ B = and [D 2 Z B ] is negative semi-definite, so 

dZ B f t 



m ij (cp'Dd) 



dt 2M 



D ylyJ Z B + —DidDjd + f'Dijd 



ft 
2M 



= m W yV Z B + W[YT ^ + f D vl d - — , 

where we have used that \Dd\ = 1 and DidD^d = 0. 
As in Lemma 14.17 of [16| . 

n-l 



tiace[D 2 d(y)] = J^y 



Kid' 

i=i 

where ft, are the principal curvatures of d£l at x, the closest point on to y. If d < R, 
then Kid < 1 and 

n—l n—1 



V — ^— < -V Ki < o, 

^ 1 - Kid ~ 4^ 



i=l i=l 

the last inequality resulting from the mean-convexity of dfl. Then 

3Zb ij n 7 1 ^" i / t-\ i ft , n 

It follows that Z B < for i < X", and so for all x e dQ and y e -B R (x) 



- tt(ai,i) < 2M<^ 



/ dist(y,an) t \ / Jx-yj t \ 



This gives us an estimate on the boundary. We complete our proof by using the 
same Z as before: 



Z(x, y, t) = u(y, t) - u(x, t) - 2Mcp ( 



x — y\ t 



\ 2M 'AM 2 

Once again, Z(x,y, 0) < 0, and when both x and y are on the boundary, Z(x,y,t) 
-2Mip < 0. If x is a boundary point and y is an interior point (or vice-versa), then 

Z(x, y, t) = u(y, t) - u(x, t) - 2M<p (^^, j^j < 
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bv grrn . 

Finally, if (x,y) is a maximum of Z at some time t < T', where both x and y are 

dZ 

interior, then as in Section EH-^- < at this point, and so Z < for all t < T . The 
estimate follows. □ 

Remark: We can use these methods to find gradient estimates for equations of more 
general form. 

For the Dirichlet problem with conditions on a ij given in Theorem 15.11 and u = 
on d£l, we can find estimates of the type in Theorem l5~4l for convex fl 

If a ij has the form Q5.5) . then we can find estimates of this type on domains that 
are merely mean-convex. 



Chapter 6 



Application of gradient estimates 
to the Neumann problem 



In this chapter we use the gradient estimate derived previously to establish the exis- 
tence of solutions to the mean curvature flow equation with Neumann boundary con- 
ditions 



where O c R n is a compact, open convex domain with C 2+a boundary <90, and u e 
C(H). The outward unit normal on the boundary is v. 

This extends Huisken's result in [TJ) showing the existence of smooth solutions to 
(fQ) for initial data with greater regularity. 

Theorem 6.1 (Huisken). Let n be a bounded domain in R n with dn e C 2+a . If 
u e C 2+a (Q) satisfies D u u = on dQ, then \6.1) . 16.2) has a smooth solution on 

n x (o,t). 

Note that while this theorem makes no restriction on the convexity of 0, the main 
result of this chapter does. 

Theorem 6.2. Let ft c R n be a smoothly bounded, open, convex domain, and let 
u e C(TT). Then the Neumann problem \6A) has a smooth solution fort > 0, which 
converges uniformly to u ast — >• 0, at a rate dependent on the modulus of continuity 

OfUQ. 

6.1 Some remarks about changes of coordinates that 
straighten boundaries 

A similar discussion of boundary curvatures and the distance function may be 
found in Appendix 14.6 of (T6). Consider a bounded domain Q c K n with bound- 
ary dn. The boundary is said to be C k+a if for each boundary point x we can find a 




(6.1) 



D„u(x, t) = for x e dn and t e (0, T] 



(6.2) 
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(jk+a ma ppj n g f ■ ^n-i ^ wn j c h has the boundary in a neighbourhood of x as its 
graph. 
Set 

i? := i sup { r : If dist(x, 00) < r then x has a unique closest point x e <90 } . (6.3) 
If is convex, then we can take 

R := - inf — — , 

l<i<n-l v ' 

where are the principal curvatures of <90. 

On balls B R (x ) centred on the boundary, we introduce a change of coordinates 
* : Br — > 50 x [-.R, .R] such that if the new coordinates are denoted y = (y, y n ) = 
(y 1 , . . . , y n_1 , y 71 ), ^(x) = (y, y n ), and I : «90 — > M. n is the immersion of the boundary 
into M n , then d(x, 90) = d(x,I(y)) = y n . In other words, i(y) is the closest point to x 
on 00 and y n is the signed distance between x and y, being positive if x O, zero if 
x g 00, and negative otherwise. 

The inverse transformation is easier to work with, being given by * _1 (y, y n ) = 
/(y) + v(y)y n , where z/ is the outward-pointing unit normal to <90. 

As the boundary in the new coordinates is simply y n = 0, this is referred to as a 
boundary-straightening transformation. 




If the graph I(y) = (y 1 , . . . , y n 1 , /(y)) is a local immersion of the boundary, then 
the outward unit normal is given by 



i 



Of 



of 



and on ty(B R (x )) we have 



Oy* 

[4 + V n w fori,i = l,...,n-l 
§t + y n %^ forj = nandi = l,...,n-l; 



(6.4) 



for 



i = n, 
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so eigenvalues for [D^- 1 ] are 1 - y n Ki(l + [Df] 2 )- 1 / 2 , . . . , 1 - y n K n ^(l + (D/l 2 )- 1 / 2 , 
and lie between \ and § on B R . The curvatures of the boundary, are given by the 
eigenvalues of [D 2 f]. 

Also, second derivatives are 

' y n WW for i,j,k^n, 

^h + y n S£z for j = nandi,fc/n, 



QyiQyk tl dy l dy 
dvj 
dyi 

dy k 

for i = n,k = n. 



for i / n and fc = n, 
for i = n and fe / re, 



The smoothness of this change of coordinates is dependent on the smoothness of 
the boundary: if dn is c k+a+2 then D^ 1 is C k+a . 

When u is defined in the old coordinates on M. n , in the new coordinates we can 
define a new function 

v(y,t) -u^iy)^). 



First derivatives are related by Div(y,t) = [D^f l ] k D k u(^ 1 {y),t), and second 
derivatives by [D 2 v]ij = [D^~ 1 ] 7 p[D'^~ 1 \[D 2 u] 7n i + (D m u)[D 2 * -1 ]^. 

Putting this all together, we notice that if u satisfies \QA \ then v satisfies 

v t (y,t) = m?t([D*]Dv) ([D^D^D^ + {D s v)[DH]\ 

= a kl (Dv,y)[D 2 v] lk + b(Dv,y) for |/6f(Q fi ), (6.5) 

!>„«(!/,*) = aty n = 0. (6.6) 

where the mean curvature operator is abbreviated as m ij (p) = 5 ij - ^ifff 2 ^ , and we 
write 



a ifc (p,y): = m^([^b)[^. [jD ^ 
b(p,l/) —m^aD^Cp.)^ 2 *]^. 



(6.7) 



Once we have straightened out the boundary, we will find it useful later on to define 
a reflection p in the boundary that extends v outside fi: 

v(y,t) :=v(p(y),t) for y G Br 

where p(y) = {y\... ,y n ~\-\y n \). 

Let Q R be the intersection of a parabolic cylinder with the domain of interest: 

Qr{x , t ) = { (x, t) G n x [0, T] : x G B R (x ),t G (t - # 2 , i ) }•) 
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When v satisfies <6.5| on Q R , v will satisfy 

v t (y,t) = v t (p(y),t) 

= a kl (Dv(p(y),t),p(y))[D 2 v(p(y), t)] kl + b(Dv(p(y),t), p(y)) 

= a kl ([Dp]Dv(y,t),p(y))[Dp]l[Dp]i[D 2 v(y,t)] l , + b([Dp]Dv(y,t), p(y)) 

= ~a kl (Dv,y)[D 2 v] kl + b(Dv,y) (6.8) 

on B R x (t! - R 2 , h) where [Dp] = diag(l, . . . , 1, -y n /\y n \). 

The regularity of the coefficients of the reflected equation is estimated: 

Lemma 6.3. Ifv is a C 1+a function on Q R , with D n v(y, 0) = 0, then the coefficients for 
the reflected equation 

&{Dv,y) = a kl ([Dp]Dv,p(y))[Dp]l[Dp]l 
b(Dv,y) = b([Dp]Dv,p(y)) 

satisfy Holder estimates 

\a l3 \ a -B R x{tx-B?) < 2 \ al3 \&;Q R (! + \Dv\atQ R ) i 
Ha;B«x(ti-fl2) < 2|6|ci;Q H (1 + I^I^Qh) 

for some < a < 1. 



Proof: In general, if a function is defined piecewise on a convex domain £7 divided 
into Ui and C/ 2 = U\U U 

h(x)= \hi(x) zeUi 
\h 2 (x) x£U 2 

and is continuous across any shared boundary Ui n U 2 then if /i x is C a on J7i and h 2 
is C a on C/ 2 , it follows that h is C Q on U = U x u J7 2 - 

We can see this by letting xi and x 2 be in U\ and ?7 2 respectively. We can find 
a point z in the shared boundary U\ n F 2 directly between the two, with |y - x\ = 

\y — z\ + \z — x\. 
Then 

\h{x)-h(y)\<\h{x)-h{z)\+\h(z)-h{y)\ 

= \h l (x)-h 1 {z)\ + \h 2 {z)-h 2 {y)\ 
< C\x - z\ a + C\z - y\ a 
= C (s a \x - y\ a + (1 - s) a \y - x\ a ) 
= 2C\x-y\ a 

for s = \x — z\/\x — y\ < 1. 

This observation applies to both v and Dv — as D n v = on the boundary, Dv = 
[Dp][Dv] is continuous across the boundary, even though [Dp] itself is not — and so 
Dv is C a . 

It is clear that b(Dv, y) = b ([Dp]Dv, p(y)) = b (Dv, p(y)) is continuous, and b shares 
the same regularity as b. 
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As 



a ij (Dv,y) 



a ij ([Dp]Dv(y,t), p(y)) for 1 < i, j < n - 1 or i = j = n, 

-$rja ij {[Dp]Dv(y,t),p(y)) , for i ^ n and j = n or vice-versa, 



we only need to check whether the terms in the off-diagonal block d in are continuous. 
These are given by 



~a in (Dv,y) 



' n a m ([Dp]Dv(y,t),p(y)) 



\y n 



\y 



y -m kl ([D*][Dp]Dv) [D*]i[D*]f 



-y n 



t _ g^f ([DnDp]Dv) a ([DnPp]Dv) p [D*]i[D*]? 

\y n \ \} k k l + \[D^][Dp]Dv\ 2 

yn_ ^^^[D^llDp^D^D^DpY^DsvjD^D^ ^ 
\y n \ V l + \[D*][Dp]Dv\ 2 J 

n i / n ~i n—1 

\y-\l + \[DnDp]Dv\\^ p 7 J ^ 



^ J2[DWaD 7 v[D^D n v - jL- Y^[D^D n v[D^D p 

+ [D^D n v[D^D n v] [D*]%[D*]Y, 



\y \ — V' , 

|y 1 7=1 |y 1 ^=1 



(here there is no summation over n). 



Between the third and the fourth line, we have used that [D^^ 1 ^ is tangent to the 
boundary while [L>^ _1 ] n is normal to the boundary (see equation i|6.4| ). so for i / n, 
we have EaP* -1 ]*^* -1 ]^ = . It follows that [£>tf]£[D¥]g = 0. 

In the last step, the second, third and fourth terms are zero (and so continuous) on 
the boundary, as D n v(y, 0) = 0. The first term is zero due to the presence of 

£(P*]fiV) W = o, 

since p^n. 



So, both a ij and b are continuous across the boundary. 



It follows from the first observation that the Holder constant of a ij on B R (y{) x (t\ - 
R 2 ,h) is the same as that of a ij on QR(yi,h); and if we consider a ij (Dv(y),y) as a 
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function of y, then we find that 

\a* i (Dv(-),-)\ a <2\a i *(Dv(-),-)\ a 

_ 9 \a ij (Dv(zi),z 1 )-a i ^Dv(z 2 ),z 2 )\ 
— — sup . . 

zi,z 2 \Z\ ~ Z 2 \ a 



1 

< 2 sup i 1 

21,22 \ z l — z 2\ 



\a»(Dv(z 1 ),z 1 )-a^(Dv(z 1 ),z 2 )\ 

+ \a i ^Dv(z 1 ),z 2 )-a^(Dv(z 2 ),z 2 )\ 



< 2 S up\a lJ (Dv( Zl ),-)\ a 

Z\ 

+ sup \a lj (-, z 2 )\ c i sup | - — —\Dv(zi) - Dv(z 2 )\ 

Z2 Z\,Z2 \ Z 1 Z 2\ 

< 2\a ij \ c i (1 + |£H;| a ) . 

Similarly, if we consider b(Dv, y) as a function of y, we find that 



\b{Dv{-) 1 -)\ a <2\b{Dv{-),-)\ 

< 2|6| c i (l + \Dv 



a 



□ 



6.2 Existence of solutions with continuous initial data and 
Neumann boundary conditions 



We begin our proof of Theorem E2] by approximating the continuous initial data by 
mollified functions that will satisfy the requirements of Theorem l6.1l being smooth and 
satisfying the Neumann boundary condition. 

Lemma 6.4. There exists an approximating sequence u e C°°(Tl) with u e ->■ u in 

C(tt), D u uq = on dtt, andu ei > u t2 > u whenever ei > e 2 . 

Proof: Let B R be a ball centred on the boundary. We work in the new coordinates 

on V(B R ), and write v (y) = u (1'~ 1 (y)). 

Remembering that D denotes the extension by reflection of v Q , define the mollified 
function 



v o(y) '■= Ve * vo (y) = / Ve(y - z)v {z)dz, 

Jza^{B R ) 

where we use the usual mollifier 



otherwise, 



and7? £ (z) = ir?(f). 

This approximation has all the usual qualities of mollifications: vq g C°°(^(B R ) e ), 
where *(B R y = {y e V(B R ) : dlst(y,d9(B R )) > e} ; and since v e C(V(B R )), 
-> v uniformly on compact subsets of V(B R ). 
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In addition, each satisfies the Neumann condition D n v^(y) = for y e 90 n 
*(B R y, since 



D nVo(y) = D n r) e (z)v (y - z)dz 

JzeB e (0) 

= / 6 B e (0) *W D MV ~ Z ^ Z + / eBe (0) ^)A^o(</ " 



2™>0 ^"<0 

Recalling the relationship between the reflected and original functions, v (y - z) = 

voiy 1 - z 1 , . . . , - \y n - z n \), we observe that when y e 90, y n = and D n v {y - z) = 
^Dnvoiy 1 -z 1 ,... , -\z n \). Consequently, 

D n v e (y) = /_ , n , Ve(z)D n v (y 1 - z\ . . . , -\z n \)dz 

z n >0 

- J z ^ {Q) Tit(z)D n Vo{y x -z\...,-\z n \)dz = 

z n <0 

as the mollifier has the symmetry rj^z 1 , ... ,z n ) = ^(z 1 , ... , -z 11 ). 

This is only a local approximation, but in the next step we extend it to the entire 
domain, taking care to preserve the Neumann boundary condition. 

Let the set of boundary-centred balls { B R (xi) }»=i ) jv be a finite cover of the bound- 
ary 90 with the property that the set of balls of half the radius { B R / 2 {xi) }i=i,jv is also 
a cover. On each ball B R (xi) we can define the approximation vq { := Vq as described 
above. 

Now, define a new cover of by the sets 

W i := {x E B R ( Xi ) : if *(x) = (y,y n ), then y E V(B R/2 ( Xi )) and \y n \ < | }. 

The cover is completed by W° := { x e : dist(x, 90) > i2/4 }. 

Note that B R/2 (xi) c c B R {xi) and so this is indeed a cover; also, is 
defined on W\ On we define the usual mollification with no reflection, which we 
call ul^ 

Let ^ be a partition of unity with respect to the sets {B R / 2 {xi) n 90} which cover 
the boundary; that is, < ^(x) < 1 for x e B R/2 (xi) n 90, & e C^(B r/2 n 90) (that 
is, compactly supported with respect to 90), and = 1 f° r a " x G 

In the new coordinates on B R (xi), we could write ^ = ^(y 1 , . . . , y"" 1 ), since ^ is 
defined only on the boundary. We can extend ^ to all of B R / 2 (xi) by setting £;(x) := 

i i ^(x)\...,^(xr~ i ). 

Let C : K ->• K be a smooth cut-off function satisfying 



C(«0 



1 ldl<# 

|d|>|. 



We will set ((x) ■- C(d(x,90)) where d(-,90) is the signed distance function. 

Now, we claim that the functions t,i(x)((x) and 1 - ((x) are a partition of unity with 
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respect to the sets {W { n 0} and W°. Firstly, all functions are smooth and compactly 
supported on their respective domains (but they are not zero on the external boundary 
dfl). Secondly, if x e 0, then 

^&(x)C(x) + (l-C(z)) = l- 

i 

This is because if dist(x,9fi) > R/2 then ((x) = 0, while if dist(x,dU) < R/2 then 
x has a unique closest boundary point x . In the latter case ^(x) = £i(x ) and 

EMx)((x) + (l - CO*)) = C(x)EM^) + (l - C0«0) = l as & is a partition of 
unity on the boundary. 

This construction ensures that D u (((x)£i(x)) = when x e <9S1 and so if we define 
our global approximation as 

N 

u e (x) :=^i(x)C(x)v e 0ti ^f(x)) + (1 - C(x)H fi (x), 

i=l 

we find that u e Q — >• u uniformly in (7(0), and each uq e C°°(f2) satisfies -D^uq = on 

We can ensure that this sequence is monotone in e, in the sense that u € q(x) < 
u e 2 (x) whenever e 1 < e 2 by restricting to a subsequence and off-setting if necessary. 
□ 

The result of Huisken mentioned at the start of this chapter now implies that there 
is a smooth solution u e to (16.1 > with u e (-,0) = u e Q . 

Lemma 6.5. The approximate solutions u e have a uniform height bound 

sup \u e \ < sup |«o|- 

f2x[o,T] n 

Proof: As the mollification u e is created by a local averaging of no, 

sup \u e (-, 0)| <sup|«o|- 
n n 

Suppose at some time t > and point xi, u e equals \u \. From the Comparison 
Principle (Theorem I2.2}. x\ can be assumed to be a boundary point. The Neumann 
condition D u u e = implies that Du e = and so [D 2 u e ] is negative semi-definite at this 

point; it follows that < and so u e is not increasing at this point. □ 
This height estimate is of course also an oscillation bound 

\u e (x,t) - u e (z,t)\ < 2\u \. 

We are now in a position to use the gradient estimate of Theorem Ol For some 
T > 0, there is a gradient bound 

\Du € (x, t)\ < L(t) for t e (0, T) and x e n, (6.9) 

where L(t) and T are dependent on n and sup n \u \. 
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Lemma 6.6. Higher derivatives ofu e are uniformly bounded on the interior, with 

/ \ — 

\D k u%r x{tl} <c(n,k,L(t ))(^ + j—j-j 2 (6.10) 

for ti > t > and all k = 1,2, ... , where Q r is the interior set{x££l: dist(x, dO,) > 
r }. 

Proof: This is an application of the Ecker-Huisken interior curvature estimate de- 
scribed in Theorem l2?6l originally in the paper (T4| . 

We apply it to the interior of Q (with 6 = and k = m + 2) to find bounds on all 
higher derivatives. □ 

This estimate provides no information as we approach the boundary. However, 
our uniform gradient bound L(t ) ensures that the evolution equation is uniformly 
parabolic, since for t > t , 

As we have uniform parabolicity for strictly positive times, extending regularity up 
to the boundary is a routine application of known results. This is the subject of Lemma 
1671 - LemmaEH 

We begin by showing that a function with a Holder estimate on the boundary of a 
region, and a strictly interior gradient estimate, has a global Holder estimate. We plan 
to apply this to finding a Holder estimate for the gradient Du. 

Lemma 6.7. Let £1 be a convex domain. ff/:fix[0,T]->M has a Holder oscillation 
bound on the boundary 

osc / < C(t )r a for all x Q £ 80, andt > r 2 , 

Qr(xo,to) 

where C(t) is non-increasing in t; and gradient bounds on the interior 

1/2 



and 



' / '^ t "^( d.st(x!^ + I 
then we can find an a' > such that for all x, y e fi ands, t e (0, T], 

\f( x ,t)-f(y,s)\<C(\x-y\ + \t-s\ 1 /y 

where C depends on min{t, s], diam Q, c and a, and a' depends on a. 
Proof: We split the difference in the obvious way 

\f(x,t) - f(y,s)\ < \f(x,t) - f(y,t)\ + \f(y,t) - f(y,s)\, (6.1 1) 
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and look at the first term. 



Without loss of generality, set d = dist(y, dQ) < dist(a;, dtt), and y to be the 
closest point to y in 89, so that \y - y \ = d. 



If we are close to the boundary, so that d < \x-y\, then 

\f(x, t) - f(y, t)\ < \f(x, t) - f(y , t)\ + \f(y , t) - f(y, t)\ 

< osc / + osc / 

Q\ x -y \(yo,t+\ x -yo\ 2 ) Q\ y -y \(yo,t+\y-yo\ 2 ) 

< C (t +\x- y Q \ 2 )\x - y \ a + C(t +\y- y \ 2 )\y - y \ a 
<C(t) [\x-y\ a + \y-y n+C(t)\y-y \ a 
<2C(t)\x-y\ a + C(t)\y-x\ a 

<ci\x-y\ a 

where we have used that \y - y \ = dist(y, 89) < \x - y\, and have set ci = 3C(t). 



If we are further from the boundary, so that d > \x — y\, then for some e e (0, 1), 

\f(x,t) - f(y,t)\ = \f(x,t) - f(y,t)\< \f(x,t) - f{y,t)\ 1 ^ 

f 



< 



osc 

'\x- yn \(yo,t+\ x -yo\ 2 ) ' 



+ 7 ) \ x ~v\ 



1-e 



&ist{y,d9) 2 t 

< [C(t+\x-y \ 2 )\x-y Q rYc 1 ^(^ i + -^ " \x - y\ 1 ^ 

< C{tf [\x - y\ a + \y- y \ a Y c 1 ^ + t^' 2 ) \x - y\ l ~* 

< 2 e C(t) e (d ae+e - 1 + d^fc-V'A \x - y\^ 

< c 2 \x - y\^ e , 

setting e = l/(a + l) SO that ae + e - 1 = 0, and c 2 = 2 e C(t) e (l + (diam9) at t^ 1 ^ 2 ). 



Now consider the second term of (|6.1 1f . and suppose without loss of generality 
that s < t. As before, set d = dist(y, 89) and y e 89. 



If we are close to the boundary, so that d < sji - s, then 

\f(y,t)-f(y,s)\< osc / 
< C(t)|t-s| a / 2 
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Otherwise, if d > y/t — s, then 



\f(y, t) - f(y, s)\ = \f(y, t) - f(y, s)\" \f(y, t) - f(y, s)\ 



< 



osc / 

}d(yo,t) 



1 1 

d 2 " + ~ S 



1-fi 



1-1" 



< [C(t)d 

< c{tyc l ~» + \ \ t - si 1 



< cAt 



|(2-2yu)/2 



where /i = 2/(a+2) so that a//-2(l-/i) = 0, and c 4 = C^c 1- ^ (l + (diam ft)"^" 1 ) . 

We find the final estimate by choosing C = sup{ci, c 2 ,c 3 ,c 4 } and a' = min{l - e, 
2 - 2^} = min{a/(a + 1), 2a/ (a + 2)} = a/(a + 1). □ 

Lemma 6.8. The gradient Du e is Holder continuous, with bound 

\Du e \ ajCt/2 < C 



onVtx [t, T], fort > t > and some a > 0, where C = C(n,L(t ), \t — t \,Q, \8Q\ 



Proof: We can use the Holder gradient estimate near a flat boundary from Theorem 
IA.1I but we will need to work locally with v £ (y,t) = « e (\I/ _1 (y),i) in the flat-boundary 
coordinates. 

The gradient bound Q6.9| for u e implies that 



\Dv e \ Qr < L(t )\D*- l \ Br < L(t )\dn\ 



c 2 



on the cylinder Q r (yi,h) for some r < R to be chosen later, and where yi e 89, and 

ti > t + r 2 . 

On Q r , v € satisfies the evolution equation \6.5\ . To check that the coefficients 
and b (given by 1/5.7)) satisfy the conditions of Theorem I7Q1 we note that: a lk is 
uniformly parabolic, since for £ e W 1 , 



> 



1 



> 



> 



1 + \[Dt!]Dv e \ 2 

1 A 2 

1 + A 2 - ^ 



lei 2 



D iei 2 



l + 4L(t ) 2 \3 
> c(L(t ))|e| 2 

where X[ D ^ and A [D ^ are the smallest and largest eigenvalues of [DV], which are 
bounded between 2/3 and 2 on B R ; \a lh \ is bounded above, as 



\m ij ([D^}Dv e )[D^] l AD^]f\ < |A 
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and a lk has bounded derivative with respect to the gradient, for if we write q = [D^]p, 
then a lk (p,y) = m'J (?)[£>*]} and 

daPHp, y) r . , k drrii{g) dgP 



i + M 2 1 (i + kl 2 ) 2 



< 4| J D^| 3 

< 4|<9ft|^ 2 . 

It is straightforward to bound the lower-order term in the equation — 

\m ij {[Dm}Dv e ){D 2 m} s ij D s v t \ < L(t )\D 2 ^\ < L(t )\dn\ c3 . 
Finally, we need an oscillation bound smaller than a for a lk (p, •) on Q r , but since 

da lk (p,y) am«([D*]p)5([D*]p) , r _. Tl , r _. Tlfc w ^ / rriTl , mTlfc \ 

Q y = — g q i dy [ m U m \ k i+ m3 Q- y ([mwmt) 

< C\D 2 ^\\D^\ 

<c\an\ 2 c3 , 

we find that oscg r ( yi)tl ) a iJ '(p, •) < rC|Z? 2 ^| By choosing r small enough, we can 
ensure that this is less than a. 

Now Theorem lATI implies that for all s < r there is some a' > so that 

oscDv 6 <c(-Y (oscDv 6 + L(t )\dfl\ C 2r 

Q s \T / \ Qr 

<c^Y {L(t )+L(t )\DH\r) 
< cs a 'r- a '(l + r). 

where c(L(t ),n, \D% \D 2 V\). 

This boundary oscillation estimate for Dv e on Q s {yiM) for a\\ s < r < R and 
t\ > t + f 2 , together with the interior gradient bounds given by Lemma l6~6l means we 
can use Lemma l6~71 to give a global Holder bound for Dv c and hence for Du e . □ 

Lemma 6.9. We can find bounds foru e in H 2+a x (*, T]) 
fort > t > 0, where C = C(n,L(t ), \t - t Q \, \dCl\ C 2, a). 

Proof: To establish this, it is possible to use boundary estimates for the Neumann 
problem, but instead our approach is to use the reflection v e on B R (y\) x [t ,T] — a 
domain that extends beyond £1 — which satisfies the reflected evolution equation (|6.8| . 
and apply the interior estimate from Theorem fA4l 

We need to check that equation \6.8) satisfies the conditions of Theorem fA4l 
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In Lemma 1631 we showed that the coefficients in equation d6.8fr have regularity 
estimates 

\a ij \a;B r x(ti-r2) < ^\a ij \chQr ( l + \ Dv %;Q r ) 

\Mcr,BrX(b.-r*) ^ 2 \ h \c\Q r i 1 + \ D V%;Q r ) 5 

and in Lemma|6i8]we found a uniform global bound for \Dv e \ a ^/ 2 . nx ^ iT ) for t > *o- 
Our gradient estimate ensures that a is uniformly parabolic for t > to- 
Applying Theorem |7V4] results in the bound 

\v e \2+ a ,l+a/2;Q R/3 < C \v%Q 2R/3 , 

where C is dependent on the dimension n, R (which is determined by \d£l\ C 2), the 
ellipticity constant of a ij (dependent on L(t )), the Holder exponent a, and the bound 
on the a norm of the coefficients (which is bounded by \Dv e \ aa / 2 ). That is, C = 
C(n,L(t ),\dn\ C 2,a). 

We can repeat this over the entire boundary, and together with the interior estimate 
i6.10f . this gives us the claim. □ 

Lemma 6.10. The sequence of approximate solutions u e converges 

Hm\u e - w|2+ a ',l+a'/2 -> 

to some u £ H 2+a < on tlx (t , T), for all t > 0. 

Proof: The uniform H 2+a (Q x (t , T)) bounds on the u e ensure that there is a conver- 
gent subsequence (for a slightly smaller a' < a); the disjointness of initial data u e (and 
hence the disjointness of u e (-,t)) implies that the entire sequence must converge, and 
so this limit is unique. The limit u is in H 2+a > (ft x (to, T)) (and is C°° on the interior, by 
virtue of the interior estimate in Lemmal6!6|. 

It also satisfies D v u = on the boundary, and so is a solution to the Neumann 
problem given by Q6.1|.D 

Lemma 6.1 1 . This solution u converges to u in C(O) as t — >• 0. The convergence is 
uniform in time, and ifu e C a (J7) thenu e H a (fl x [0, T]). 

Proof: Let 5 > be fixed. Our aim is to show that we can find t 5 so that for all 

t G (0,t S ), 

sup\u(x,t) — uq(x)\ < 5. 

The modulus of continuity for u is 

u(r) := sup \u (x) - u (y)\. 
\x-y\=r 

Recall the approximate solutions u 6 , converging uniformly 

\u e (-,t) - u(-,t)\ -> as e 
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for all t > 0. The approximations at the initial time have (at least) the same modulus of 
continuity as u : 

\u e (x,0) -u e (y,0)\ <to{\x-y\). 




Figure 6.1 : The modulus of continuity bounds u 



Fix z to be any point in the interior of Q. We can define a new solution to QQ) by 
off-setting u e around u (z): 

w e (x, t) := u e (x, t) — u e (z, 0) + uo(z), 

SO that w e (z,0) = u (z). 
Then, for t > 0, we have 

\u(z, t) — uo(z)\ = lim \u e (z, t) — uo(z)\ 

= lim \w e (z,t) +u e (z,0) - 2u (z)\ 

< lim (\w e (z, t) -u (z)\ + \u e (z,0) - u (z)\) . (6.12) 

To estimate the first term of this, we observe that every w € (-, 0) is inside the 'enve- 
lope' given by the continuity condition, 

— L0(\z — x\) + Uq(z) < W e (x, 0) < U>(\z — x\) + Uq(z). 



Above and below this envelope, we can place two spheres of radius r centred at 
(z,u (z) ± [r + uj(r)]). At t = 0, the spheres and the graph of w e (-,0) are completely 
disjoint. The spheres are also completely disjoint from the graph of u . 

The evolution of spheres under mean curvature flow is well-known — the centre 
remains fixed and the radius shrinks from the initial radius r(0) = r , with 



r(t) = \Jrl~ 2nt 
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Figure 6.2: Spheres above and below w e or u 



until the sphere disappears at time t = rl/2n. 

The parts of these spheres closest to the graph of w e — the lower part of the upper 
sphere and the upper part of the lower sphere — are 

S + (x,t) := u (z) +r +w(r ) - yV^) 2 — \x — z\ 2 , 
S~(x,t) := u (z) - r - u;(r ) + \f r(t) 2 - \x - z\ 2 . 

Suppose that one of these spheres and w e first touch at some time t > 0. From 
the Comparison Theorem l2~2l we know that at this time there must be an intersection 
occurring on the boundary of say at x\ e dVt (this doesn't rule out other intersections 
occurring simultaneously on the interior). 

This intersection on the boundary is an extreme point of S + /~ (•, t) - w e (-,t) (either 
a minimum of S + - w e or a maximum of S~ - w e ). 

Therefore the sign on the outward derivative of the intersecting sphere at this point 
is known — either 

D v (5 + (xi, t) - w'ixi^)) < and so D v S + (x x ,t) < 0, 

or else 

D u (S~(xi,t) - w e (x 1 ,t)) > and so D v S~(xi,t) > 0, 

where we have used that w e satisfies the Neumann condition D v w e = on the bound- 
ary. 

On the other hand, we can explicitly calculate the gradients of the spheres — 

D>S+(x,t) = 2 f~ Zl) D iS -(x,t) = ; 2 } Xl ~ Zl \ 2 . 

y?r z — \x — z\ A \Jr z — \x — z\ z 

The convexity of ft means that for any z e Hand x\ on the boundary, v{x\—z) > 0, 
and the inequality is strict if z is in the interior of £1. Hence the sign of the normal 
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derivative is known — 

D u S + (x 1 ,t) >0, D v S~(x u t) < 0. 

Between these two observations, it must be the case that the intersecting sphere 
has a flat normal gradient — 

D v S + (x 1 ,t) = or D v S~(x!,t) = 0, 

and so v-{x\—z) = 0, which in turn implies that z is on the boundary of 0, contradicting 
our original assumption that z was an interior point. It follows that such spheres, 
centred on interior points, never touch the graph of w e for the duration of the spheres' 
existence, until t = r 2 /2n. 

In particular, above the point z the surfaces move by no more than r + u>(r) in the 
time t e (0,r 2 /2n). We can choose r > so that 6 = r + w(r), and a corresponding 
t s = r 2 /2n, ensuring that 

\w e (z,t) - u Q {z)\ <r + f(r) = 5ior0<t<t 5 , 

where t is dependent on 5 and u alone. 

This estimate is independent of z and e, so 

lim \w e (z, t) — uo(z)\ < 5. 

We still need to estimate the second part of equation i|6.12| . \u e (z,0) - u Q (z)\. 
However the convergence here is uniform (for z e ft), so that 

sup \u(z, t) — uo(z)\ < 5 + sup lim \u e (z, 0) — uo(z)\ = 5 

for < t < t 5 , and so u is in C([0, T); C(Q)), with «(-, 0) = u . 

This is also an estimate for the smoothness in time of the convergence; we can 
consider 5 = S(t), with 

sup \u(z, t) — Uq(z)\ < S(t) 
zed 

by setting t = r 2 /2n, so that r = V2rd and thus 5(t) = y/2nt + u{\/2ni). The con- 
vergence in time that this gives is at best like t 1 / 2 — which is in concordance with the 
result for initial data with a Lipschitz bound. In that case, u(r) = |r| and the conver- 
gence is C 0+1 / 2 in time (Theorem 3.5 of l20l V 

In the case that the initial data has a Holder gradient bound, u(r) = \r\ a , then the 
convergence to the initial data is as t 1 / 2 + W 2 ~ t a l 2 . □ 

Remark: While we have mined the rich theory arising from mean curvature flow 
to find this result, there are similar results for other equations of the type studied in 
Chapter|Sl and we expect to be able to find similar existence results. 

In particular, one can find short-time existence results for anisotropic mean cur- 
vature flow with a zero Neumann boundary condition and continuous initial data on 
Q x [0, T], for convex £1 
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Existence of solutions to the 
Dirichlet problem for mean 
curvature flow 



In this chapter we use the gradient estimate to establish existence of solutions to the 
Cauchy-Dirichlet problem with zero boundary data. 

Theorem 7.1 (Existence of solutions to the Dirichlet problem). Let be a domain 
in M. n , with C 2 boundary 90 that has non-negative mean curvature. Ifu e C°(U) and 
u = on <90, then the problem 



has a smooth solution fort > which converges uniformly to u as t — > 0. 

The existence of solutions to the mean curvature flow problem with prescribed 
boundary values was considered by Lieberman in |24) (and by Huisken in (TJ), where 
the long-time behaviour of solutions was also studied). 

Lieberman considered time-dependent boundary data u e H i+Ct (P(£l x [0,T])), 
with a Lipschitz bound (in time) on dfi x [0, T]. The following theorem may be found 
as Theorems 12.10 and 12.18 of l25l 

Theorem 7.2 (Lieberman). Let c W 1 be a domain with C 2 boundary, and let u be a 

function defined on the boundary with u e Hi+ Q {V(Vt x [0, T])) for some < a < 1. 
Then if the mean curvature H of dQ, is non-negative, there exists a solution to i|7.1| 
with initial and boundary data 




(7.1) 



u(x, t) 

U(;0) 



forx£dQ, t>0 



u(x,t) = uo(x,t) for(x,t) e V{n x [0,T]). 



Moreover, such a solution satisfies 



[Du]f3 < c, 
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where c and depend on |«o|i+a,a/2- 

The proof of Theorem 17. 11 is very similar to that for the Neumann problem, Theo- 
rem HiU We will use the boundary-straightening change of coordinates described in 
Section l6~TI and the corresponding v (y, t) ■- no(^ _1 (y), t). 

Lemma 7.3. There exists an approximating sequence uq e C°°(U) with uq = on d£l 

andu e -> u in C{£1). 

Proof: Let B R be a boundary centred ball, where R is given by Q6.3L We will define 
a local approximation on B R and then put similar local approximations together to give 
a global one. 

In Section l6~TI we defined v to be a reflection across the boundary; this time, we 
let v be the odd reflection over the boundary 

U o{^Hy,y n )) ify n <0, 

if y n = 0, (7.2) 

-u {m- l {y,-y n )) ify»>0. 

Mollifying this in the standard way 

Voiv) ■= r]e * v = Ve{y - z)v (z)dz, (7.3) 

Jz^(B R ) 

we note that vq — ► v uniformly on subsets of ^(B R ), and we can check that if y n = 
0, then v t (y) = 0. Returning to the original coordinates, set u\{x) •— v^(x)) on 

TLr}B R (xi). 

Now, cover d£l by N such boundary-centred balls B R (xi), . . . ,B r (xn) on which 
are defined approximations u\, . . . , u e N . Complete the cover of ^ by the set £i R l 2 -.= 
{x £ £1 : dist(x,0O) > R/2}. On this interior set let u e N+1 be the usual mollification of 

wo- 
lf {£,i}i=i,N+i is a partition of unity with respect to these sets, then the sum u e = 

Ya=i & u i converges uniformly to u on U and also has u e = on <9fL 

We can restrict this to a subsequence u \ where \u - «o I < 2 ~* (b ut retain the 

notation u e for the subsequence). If we off-set each member of the subsequence, by 

replacing u ei by u ei + 3(2~ J ), then this is a completely disjoint subsequence that still 

converges to u as i — ► oo. □ 
As the boundary values 

u £ = 4 on fl x {0}, 
u e = 0on 9fix [0,T], 

are in H 1+Cc on x [0, T]), TheoremESJensures that a solution with these boundary 
values for (73} exists. Denote these approximate solutions by u e . 

The gradient estimate derived for the Dirichlet problem in Theorem l5.4l implies that 

\Du € (-,t)\ < Ci\/t(l + t)exp(C 2 /t) :=L(t), 

for constants C\ and C 2 dependent only on oscrtg < oscu . 
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Lemma 7.4. The approximate solutions u e have a Holder gradient bound 

\u e (;t)\ C ° <C, 

fort > t > 0, where C is dependent on L(t ),\h - 1 \, n, \dQ\ c n, and diamQ. 

Proof: If we revert to v e in the straightened-boundary coordinates satisfying \<o.5\ on 
Qr, Theorem lA~2l aives an oscillation bound for the gradient on the boundary-centred 
cylinder Q r {xi,h) for *i > *o + f 2 and r < R — 

osc Dv e < cr a , 

Qr 

where c and a depend on n, \a ij \ C i, \b\, A aIJ and A a a, where a ij and b are the coeffi- 
cients of <6.5| . as in <67j. These last four are in turn dependent on L(t ) and \d£l\ C 2. 
On the interior we have the bounds on higher derivatives given by Lemma l6~6l 

m+l 

\D-»u-( X ,t)\ < c(„, m ,L M ) + j^j 2 ■ (7-4) 

With m = 0, this is a gradient bound for Du. 

These can be linked together using Lemma|6i7]to find a global Holder bound 

|£> u U';Q r (xi,ii) < C, 

where a' = a/ (a + 1) and C additionally depends on \t - ti\ and diamfi. □ 

Now, if we consider the approximate solutions to begin at some time t\ > t with the 
initial data u e (-,ti), the uniform Holder gradient bound on u e (-,ti) means that Theorem 
I7.2l aives a Holder gradient bound 

\ ue \l+P;Clx[ti,T] - Cj 

for all ti > t > 0, where C and (3 depend on L(t Q ),\ti - 1 \, n, \dQ\ C 2, and diamO. 

With these uniform estimates for positive times, there must be a convergence sub- 
sequence in H 1+f3 i(0, x [h,T]) for some (3' < (3. If we off-set the initial data as men- 
tioned in Lemma l731 the convergence of the subsequence implies the convergence of 
the entire sequence to a limit u, defined on all Q for t > 0. Finally, the interior bounds 
i|7.4|i mean that on the interior, u is smooth. 

We now need to show that u(-, t) -> u as t ->■ 0. 

Lemma 7.5. As t — > 0, 

sup \u(x, t) — uq(x)\ — > 0. 
Furthermore, u has a modulus of continuity in time dependent only on that ofu . 

Proof: As in the proof of Lemma lOTI fix z to be any point in the interior of Q, and 
set 

w e (x, t) := u e (x, t) — u e (z, 0) + uq(z), 

so that w e is a solution to d) with w € (z,0) = u (z) and w e (x,t) = -u e (z,0) + u (z) 
on the boundary. 
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Let lo be a modulus of continuity for u and hence for w £ (-, 0), so that 

— Lo(\z — x\) + Uq{z) < W e (x,0) < U>(\z — x\) + Uq(z). 

Above and below these two bounds, we can place two spheres of radius r centred 
at (z, u (z) ± [r + v(r)]). At t = 0, the spheres and the graph of w e (-, 0) are completely 
disjoint. The spheres are also completely disjoint from the graph of n ( ). 

As they evolve under mean curvature, the parts of these spheres closest to the 
graph of w e — the lower part of the upper sphere and the upper part of the lower 
sphere — are 

S + (x, t) := uq(z) + r + u>(r) — \J r 2 — 2nt — \x — z\ 2 , 
S~(x, t) := uq(z) — r — uj{r) + \Jr 2 — 2nt — \x — z\ 2 . 

Initially S + - w e is positive: suppose that t e (0, r 2 /2n) is the first time that S + - w e 
decreases to zero. The comparison principle (Theoreml2^2| means that this occurs at 
some point x\ on the boundary. On the boundary, 



dt ^Jr 2 — 2nt — \x\ — z\ 2 

(since w e is constant on the boundary) and so this cannot be the first zero point; it 
follows that S + - w € > for the duration of the sphere's existence. 
The same argument shows that S~ - w e < 0. 

We can conclude that the w e move at most by r + u(r) in the time t e (0, r 2 /2n). 

This estimate is independent of z and e, so this implies that u e (-,t) — > u and 
that u e C(U x [0,T]). Furthermore, if u e C a (Q), then u e H a (ttx [0,T]) u 
C°° (fi x (0,T]). □ 
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Gradient estimates found by 
counting intersections 

In the paper [5], Angenent proved a series of results regarding the finiteness and non- 
proliferation of the zeroes of a parabolic equation in one space dimension. 

A zero of v(-,t) is simply a point x where v(x,t) = 0. A multiple zero is a point 
where both v and v x vanish. In contrast to earlier results, Angenent did not exclude 
multiple zeroes from the zero set, defining the zero set as 

z {t) = {x £ R : v(x,t) =0}. 

In the following, z(t) is often used as shorthand for the counting measure TC° (z(t)). 

These zero-counting results have been influential in many different areas, and have 
been used for geometric flows by Angenent himself, in [6], [8], [7], and [2], the last with 
Altshuler and Giga. Many others working in the area have also used these results. 

Unlike approaches that depend more explicitly on the maximum principle, this tech- 
nique seems limited to equations in one dimension. The gradient estimates found do 
not depend on the initial gradient, but do depend explicitly on the height: the smallest 
gradient estimates are found for when the height is largest. 

This work originates in an idea of Ben Andrews; also, this approach to finding 
gradient estimates has been independently used by Nagase and Tonegawa in the 
forthcoming paper |26) . 



8.1 Counting zeroes 

The estimates in this chapter rely on Theorem D of Angenent's paper: 
Theorem 8.1 (Angenent). Let v : [x , x\] x [0, T] -> R be a solution of 

vt = a(x, t)v xx + b(x, t)v x + c(x, t)v. 
such that there are no zeroes on the boundaries 

v(xi,t)^0, i = 0,l. 
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Let a, b, c satisfy 

a positive; 

a, a -1 , at, a x , and a xx bounded; 

b, b t and b x bounded; 
c bounded. 

Then ifv t , v x andv xx are continuous on (xq,xi) x [0, T], 

• fort> 0, z{t) is finite 

• ifx is a multiple zero ofv at t then for all h <t <t 2 we have z(h) > z(t 2 ). 

Consider a fully nonlinear equation on a domain Q x [0, T], where is a connected 
subset of R, 

ut = F(u xx ,u x ,u,x,t), (8.1) 
where F is parabolic, by which we mean that 

d 

—F(r,p,q,x,t) > 
or 

for all (r,p,q,x,t) eR 3 x~fix [0,T]. 

Suppose that u and ip are smooth solutions of (H3) , with 

I ^| ) I ^x | > \ll xx I , I Ut I < C x 

and 

Then we can form the difference w -.= u - ip satisfying the evolution equation 

w t = u t - (ft 

= F(u xx ,u x ,u, x, t) - F((p xx , (p x , tp, x, t) 

= / —F(su xx + (l — s)<p xx ,su x + (l-s)<p x ,su+(l-s)<p,x,t)ds 
Jo ds 

f 1 d f 1 d 

= J q^ f (- ■ ■ )ds [u xx - p xx ) + J q^ f (- ■ ■ )ds (u x - ip x ) 

+ f T/ ( - )da(u - v) 

= A(x, t)w xx + B(x, t)w x + C(x, t)w, (8.2) 
where the omitted argument of the derivatives of F, denoted by (...), is always 



72 



(su xx + (1 — s)<p xx , su x + (1 — s)ip x , su + (1 — s)tp, x, t). In the last line, 
f 1 d 

A(x,t) := / — F(su xx + (1 - s)ip xx ,su x + (1 - s)ip x ,su+ (1 - s)tp,x, t)ds 
Jo or 

f 1 d 

B(x,t) := / —F(su xx + (l-s)ip xx ,su x + (l-s)tp x ,su+{l-s)(p,x,t)ds 
Jo °P 



and 

C(x,t) :-- 



f 1 d 

I —F(su xx + (1 - s)cp xx , su x + (1 - s)ip x , su + (1 - s)</?, x, 

JO °1 



In order to use Angenent's theorem, we need to establish that: 

• A, A' 1 , A t , A x , and A xx axe bounded, 

• B, B t and B x are bounded 

• and C is bounded onflx [0, T]. 

Let/C = { (r,p,q,x,t) eR 3 x!lx[0,T]: < C x }■ 

dF 

If — — is continuous, then there are positive constants A*; and for which 

or 

dF 

< Ajc < — < Ak for all (r,p,q,x,t) e /C. 



Bounds on A and A^ 1 follow from this: 



A < A 



^ is given by 



A' 1 < A^ 1 . 



d 2 F 



o 



^ r 2 (• • • ) t SU ^* + ( X ~ S )<fet] + fl^C • • ) { SU *t + _ 



and so 



1^1 < 



<9F 



9r 



(Ci + |u xxt 

<#rftax[0,T]) • 



Similarly. 



1 <9 2 F <9 2 F 

-^-(. . . ) [SU^ + (1 - sjtfea:] + q^q-(- ■ • ) [SU XX + (1 - S)^] 

<9 2 F <9 2 F 

+ TT^ _ (- • • ) [SU X + (1 - S)^] + . . )da, 



so 



I Ac I < 



drdq 
dF 



dr 



(C\ + \u XX xi 1 Pxxx\q,x[Q,T[)- 
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A xx is given by 

-i d 3 F 



o Or 3 



(...) [su xxx + (1 - s)if xxx ] 



+ 2 



3 F 
2 r0p 



(...) [su xxx + (1 - s)ip xxx ] [su xx + (1 - s)ip xx ] 



+ 2 



3 F 



d 2 rdq 



(...) [su xxx + (1 - s)ip xxx ] [su x + (1 - s)<p x ] 
3 F 

+ 2 Q2 r Q x (- ■ ■ ) [ su xxx + (1 - s)cp xxx ] + i?l 

d 2 F 

+ -foJV ■ ■ ) l su »« + (! - + #2 ds, 

where i?i and i? 2 are combinations of terms involving second and first derivatives of 

dF 

— — respectively. Consequently, 

Or 



\A XX \ < 



OF 



Or 



C 2 (/C) 



(l^asasa!) fxxx |f2x[0,T] + Cl 



+ 



<9F 



(9/' 



C!(/C) 



l^xxxxi ^Pxxxxi "ran fxxx bx[o,T]+Ci). (8.10) 



The bounds for 5, its derivatives, and C follow in a similar manner: 



2 F 2 F 
o 0p0r ^' ' ' ^ + ^ ~ S ^ ra ^ + dp 2 ^" " ^ ^ Mx * + ^ ~ 

<9 2 F <9 2 F 



so that 



l#tl < 



<9F 



Op 



(C\ + |«xa;t> <Pxxt\nx[0,T]) 



(8.11) 



while 



d 2 F 
OpOr 



2 F 



■ ■■) [SU XXX + (1 - S)(f xxx ] + ^-jO • • ) [sUxx + (1 - s)<Arx] 

<9 2 F 2 F 

+ 7nr(- • • ) + ^ " s )^-] + ttit • • • 

apag OpOx 



so that 



1^1 < 



OF 



Op 



(C\ + |« ra3: |six[0,T] + Wxxx\nx[Q,T\) ■ 



(8.12) 
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and finally 

\C(x,t)\ 



f 1 d 

/ -z-F(su xx + (1 - s)ip xx , su x + (1 - s)ip x , su + (1 - s)tp, x, t)ds 

Jo Q( i 



< 



OF 



dq 



(8.13) 



0\K 



It is clear that we will be able to apply Angenent's result to a smooth solution of a 

nonlinear parabolic equation u t = F, when F satisfies the parabolicity condition \8.6\ 

dF 

and both — — and F are C 2 on the bounded domain JC. 

or 

These conditions are not optimal — for example, in estimate <8.13| above, it is 

dF 

sufficient if — is L 1 along line segments in K — however, they are enough to allow a 

oq 

theorem for intersections of two solutions rather than zeroes of one solution. 



Theorem 8.2 (Intersection-counting theorem). Letu andup : [x ,xi] x [0,T] 

solutions of 

Ut — F{u xx , u x , u, x , t) , 
which do not intersect on the boundaries 

u(xi,t) / <p(xi,t), i = 0, 1, te[0,T]. 



be 



Ifu and areC 2 on (x ,xi) x [0,T] 



\U\ C 2, \<f\c2 < Cl, 



and if F is parabolic 



0_ 

Or 



F(r,p,q,x,t) > 



dF 

and if both F and — are C 2 on 

or 



KL={ (r,p,q,x,t) 6l 3 xfix[0,T]: M, M, |g| < ci } , 

f/7en fort > number of intersections ofu andip are finite; and if x is an intersection 
of u and cp at I then for all ti < t < t 2 , the number of intersections at t\ is strictly less 
than the number of intersections at t 2 . 



Proof: We apply Theorem (8jJ to the difference w = u - ip which satisfies equation 



8.2 Gradient estimates for equations in one space dimen- 
sion 

In this chapter we seek interior estimates for bounded solutions of parabolic equations 
on connected domains $7 x [0, T\. When two functions intersect at a single point, then 
the gradient of the one that is smaller on the left of the intersection will dominate the 
gradient of the one that is smaller on the right. 
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The main idea is that optimal regularity of u(x,t) is found by comparison to the 
solution of the same parabolic equation with initial data Co, where a is the maximal 
monotone graph 



The method can be broken into the following steps: 

• Creation of family of barriers with ip e ' s (x,t) approaching Ca(x - s) as 
t,e -> 

• Show that for all (x, t) in a subdomain of x [0, T], and for all k e [-M, M], we 
can find an s such that ip e > s (x, t) = k 

• Show that |(/' s | > M at the boundaries of fi 

• Then use the Angenent result to count the intersections of u and (p e ' s . For small 
enough e, there will be only one intersection, and so a gradient bound will follow. 

We begin by looking at a simple estimate for entire solutions on E, then find more 
specific estimates, firstly for the heat equation 



and then for a nonlinear problem. 

The following theorem says that if a solution with the step function as initial condi- 
tion exists, then it will serve as a barrier for other solutions. 

Theorem 8.3. Consider the parabolic equation 




(8.14) 



which we will refer to as the step "function". 



Figure 8.1 : The step function a 



1 



(8.15) 



u t = F(u xx ,u x ,u) 



(8.16) 
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Figure 8.2: Evolution of a step function under curve shortening flow 



where F satisfies (|8.6| . Let u be a solution of Q8.16fr on R x (0, T] fA?a? ftas a bound 

\u\ < M, 

and has a uniform gradient bound at t = 0. 

Suppose there exists a solution to Q8.16| on R x (0, T] i/i/ft/cft /'s smooth for t > 0, 
and ftas /'n/'f/'a/ condition 

ip(x,0) = (M + l)a{x) 

and boundary condition 

lim ip(x,t) = (M + l)cr(x). 

\x\—>oo 

Then there is a gradient estimate 

u x (x,t) < ip x (z,t) 
where z is chosen so that ip(z, t) = u(x, t). 



Proof: Let a family of barriers indexed by (z, r) be given by ip z > T (x, t) := tp(x-z, t+r) 
for all z e R and r > 0. Each of these satisfies Q8.16| . and is smooth on R x [0, T}. 

As n(-, 0) has a uniform gradient bound, there exists a r' > such that not only do 

n(-,0) and v3 2 ' T '(-,0) intersect only once, but also, u(-,0) and ^' r (-,0) intersect only 
once for all r e (0, r']. 

Let be fixed. For each t < r', there exists z such that ^ z,T (xi,ti) = 

u(xi,ti). 

Now, apply Angenent's theorem to w = u - (p z ' T on some region [-R,R] x [0,ti] 
containing xi and which is sufficiently large enough that for all t e [0,*i], tp z > T (R,t) > 
M and ip z ' T (-R,t) < -M. The last conditions ensure that u> has no zeroes on the 
boundary. 

As w has only one zero at t = 0, it has no more than one zero for all t; as w 
is positive at x = -R and negative at x = -R, it has exactly one zero for all t. In 
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particular the zero at (xi,t{) is the only zero, and 



forx>xi, (p z ' T (x,ti) > u(x,ti) 
for x < xi, (p z ' T (x, ti) < u(x, ti) 



from which we find a gradient estimate: 



u x (x 1 ,t 1 ) < ip z ' T x (x 1 ,t 1 ). 



This holds for all r G (0, n) and so letting r — >• gives the result. □ 

The following theorem describes an explicit barrier in the case of the heat equation. 

Theorem 8.4 (Gradient estimate for the heat equation). Let 9, = [x ,xi] and u : 

x [0, T] ->■ M be a smooth solution to the heat equation i|8.15Ji with a height bound 
\u\ < M and Lipschitz bound Lip u(-,0) < oo. 

Then for t > 0, 



This leads to an estimate for an entire solution. 

Corollary 8.5. Letu : E x [0, T] ^ M be a smooth solution of the heat equation i|8.15> 
with \u\ < M. 
Then for t > 0, 



Proof: Theorem EZJapplies on any interval [-R, R]; let R -> oo to find the result. □ 
Proof of Theorem US Without loss of generality, let O = [-d,d]. 

For some e g (0,e ), \s\ < d and N > M to be chosen later, we can define the 
barrier 



which satisfies the heat equation (|8.15| on £1 x [0,T]. As t + e -> 0, tp e > s (x,t) -> 
iVcj(x - s), where cr is the step function (|8.14f . 

Choose e small enough, such that for any |s| < d and e < e there is only one 
intersection of u(-,0) and cp e < s (-,0). This is possible as u(-,0) has a uniform gradient 
bound. 

Now, let (x,t) g O x (0, T] be fixed, and consider it := u(x,t). The bound on « 
implies that n g (-M, M). 




where 
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If we choose s = x 



^— ^Ji£rf 1 (u/N), then (f e > s (x,t) = u. Also choose 



N = M 



Srf 



\fc(d - \x\) 



n -1 



(8.17) 



With these choices, we can check that 



\s\ < \x\ + \ 1-^- l^rr 1 (u/N)\ 



x + 



't + e 



< \x\ + 



' t + e 
d — \x\ 



frf- 1 ( A frf ( ^LZM 



v^C^ ~ 1^1) 



2v / rr 



< (i; 



and that on the boundaries \ip e ' s \ > M whenever t < t , since 



\<p e ' s (±d,t)\ = NSrf (\±d-s 



t + e 



NSri 



+d-x + J — frf" 1 (u/N) 



t + e 



> N£vi \+d- x\ 



t + e 



t + e V t + e 



£ir l (u/N)\ 



> NSri (d-\x\) 



t + e 



t + e V t + e 



Srf- 1 (M/N) 



and if we use (|8.17| for (d - \x\), then this is 



l/ — Srf- 1 (M/N) - J t —£rr 1 (M/N) 



t + e 



t + e 



> M. 



We can now apply the intersection counting Theorem l8~2l with w = u — (p e ' s , with 
the previous calculation ensuring that there are no intersections on the boundary, and 
with the coefficients in equation \8.2\ given by A = 1, B = and C = 0. Since there 
is only one intersection at the initial time, there is never more than one intersection at 
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later times (in particular, for our given (x, t), there is no other intersection at time t than 
the one at x). 
It follows that 

fory>£, <p e ' s (y,t) > u(y,t), 
fory<5, f e ' s (y,t) <u(y,t), 

from which we find a gradient estimate: 

u x (x,t) < tp e ' s x (x,t). 

This holds for any smaller e > 0, so letting e -> gives the final result. □ 

This method applies to all parabolic operators for which we can find solutions that 

have the step function as the initial condition. 

When a quasilinear parabolic equation satisfies the conditions of Section l4~4l we 

can use the solutions with stepped initial data, whose existence was shown in that 

chapter. 

Let a > be in H a (IC) for all bounded KcRxlx!)x[0,T], and some a e 
(0,1). This implies that for every such K. we can find positive A^ and A& such that 
A/c < a(p,q,x,t) < Ajc, when (p,q,x,t) G K. 

Since we will be looking at bounded solutions in Q, the bound on the gradient is 
the pertinent bound on fC; we highlight this by writing: 

X(K) < a(p, q, x, t) < A(K), when \ P \ < K, (8.1 8) 

where we assume that \q\ <M,xeQ and t e [0, T\. 

Also, suppose that there are positive constants A and P such that 

a(j>, q, x, t)p 2 > A > 0, for \p\ > P. (8.1 9) 

Theorem 8.6. Letu :l!x [0, T] — >• R be a smooth solution to 

u t = a(u x ,u, x, t)u xx , (8.20) 

where a satisfies i|8.18| and i|8.19| . 
Letu be bounded, \u(x,t)\ < M. 

Let Lp s solve i8~20l on R x (0, T), with tp s {-,t) -> 2Ma(x - s) as t -> 0, where a is 
the step function (|8.14| . and where s is chosen so that u(x,t) = ip s (x, t). 

Ift < c dist(x, dnf/A (cM/dist(x, dQ)), where A(-) is given by (l8~T8l . then 

u x (x,t) < ip s x (x,t). 

That is, the gradient of u is bounded by the gradient of the barriers, at the same 
height. 

Remark 1 : We can replace a(x - s) by <x(s - x) here, in which case we find that 



u x (x,t) > tp%(x,t). 



80 



Remark 2: If a has polynomial growth, so that A(K) < c(l + K q ) for some q > 0, 
then the interior region on which we can find bounds of this form is given by t < 

c dist(x, dQ) 2+q M~ q / 2 , for some constant c. 



Remark 3: If SI = R in Theorem l8T6l then the gradient bound applies for all t e (0, T]. 

Proof of Theorem IO: 

We initially assume that 0, = [-d,d\. We will derive a gradient bound at a single 
point (0, t), and then generalize it to interior points on a general domain. 

As u is smooth there are bounds on the first derivative and on higher derivatives 

\^x\ ^ Ci, | U X x 1 1H i ^ xt i U X xx ) U X xxx j U X xt I C2- 

For < e < e << d/4, let be the standard mollification of y> 

yj e := 2Mjj £ * a, 

where a is the step function i|8.14| . 

Choose e small enough so that for all e < e , 9? e satisfies a gradient estimate from 
below: 

W\{x)\ > ci whenever < M. (8.21) 

Now, for fixed e, define a family {</? €,,s }|s|<d/2 of barriers, each of which solves \%.20\ 
on R x [0, t] (for some r to be decided later) with initial condition 

^(2,0) = v? e (x - s). 

The existence of such solutions follows from Corollary T4.20I which applies as a satis- 
fies d8~19l . 

Standard results give 



\<P e /\<c 3 (e) 

\,„e.s ,„e,s ,„e,s ,£,s ,e,s e,s i ^ / \ 



To avoid intersections of u and the barriers occurring on the boundary, we need to 
show that \ip e ' s \ > M when x e {—d, d}. 

Each barrier in the family is initially bounded above by a step function 

tp e ' s (x,0) <2Ma{x - s + e) 

and so Corollary \4. 11 I provides an estimate for x < s - e 

x 8M [At 

</' s (x,t < , —J—-2M, 

\x — s + e| V vr 

where A = A (4M /|x - s + e|). 
In particular, at x = -d, 

f a / 7 \ 8M /At 

— a — s + e V vr 
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As \s\ < d/2 and e << d/A, 



1 4 / AM \ fl6M\ 
< - and A : r < A I —r- \ 



\ — d — s + e\ d \\ — d — s + e\ J \ d J 

d 2 ir 

If we choose r = — wnere A = A(16M /d), then 

... j , 32M /AT 
^(-d,t)<_^_-2M 

< — M, 

whenever £ < r. 

A similar calculation for the other boundary point x = d gives that 

tp e ' s (d,t) > M 

when t <t. 
Let e be fixed. 

For each s e [-d/2, d/2], we can define w.= u- ip e ' s satisfying 

wt = Aw xx + Bw x + Cw 
OF 

on \—d, d] x [0,7-]. Here, A is given by setting —(r,p,q,x,t) = a(p,q,x,t) in d8.3fr . 

-B by setting — = r—a(p,q,x,t) in i|8.4t . and C by setting — = r—a(p,q,x,t) in 

oq 

E3. 



Let /C = { (p, g, x, t) : \p\ <ci + c 3 (e), \q\ < M,x £ Q,t £ [0, r] }. 

A, A" 1 , A t , A x , A xx , B , B t , B x , and C are bounded by constants dependent on a, 
c 2 , c 3 (e), c 4 (e), A (ci + c 3 (e)), A (ci + c 3 (e)) and |a| c2(/c) , as in (I8~7l-(I8.1 3>. 

Since |n| < M, and at the boundary \ip e > s \ > M, w is nei/er zero on the boundary. 

In particular, w(-d,t) = u(-d,t)-ip e < s (-d,t) > and w(d,t) = u(d,t)-(p € ' s (d,t) < 
0, so there is always at least one zero of w. The lower gradient bound i|8.21| implies 
that there is at most one zero of w at the initial time. 

Then the intersection counting theorem (Theorem l8.2^ implies there is exactly one 
zero of w for all t < r. 

In the following lemma we show that given (x, t) — or more specifically, (0, t) — we 
can find s such that (0, t) is a zero of w = u — ip e,s . We will then return to the proof of 
TheoremES 



Lemma 8.7. Let e < e and t < t be fixed. For each k e [-M,M], there exists an 

s e [-d/2, d/2] such that 

ip e ' s (0,t) = k. 

Proof: Firstly, we check that ip e ' d/2 (Q,t) < -M and <p e - d/2 (0,t) > M. Using Corol- 
lary EUJJ 

.rf/o. s 8M /At 

w £ ' d/2 (x,t) < 1 :\ 2M 

* K ' - \x - d/2 + e V vr 
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Figure 8.3: A family of barriers 



for x < d/2 - e, where A = A(AM/\x - d/2 - e|). Since \d/2 - ej" 1 < 4/d, at x = we 
have 



It can similarly be shown that (p e ~ d / 2 (0,t) > M. 

As v9 e ' s (-,0) is continuous in s, and a is a continuous operator, tp e < s (0,t) is also 
continuous in s. In particular, {<£ e ' s (0,i)}| s |< d / 2 is onto [— M, M]. □ 

Continuation of the proof of Theorem I8.6t Lett < t be given. From the previous 
lemma, there exists s such that (p e ' s (0, t) = tt(0, i). This is the only intersection point of 

u and if"-' 3 , and so 




C — M. 



fory>0, <p e ' s (y,t) > u(j/,t) 
fory<0, ¥> £,fl (y,i) < u(y,t) 



from which we find the gradient estimate: 



u x (0,t) = lim 



u(y,t) - u(0,t) 

y 



< lim 

y^O 



y 



This estimate holds for all e g (0, e ]. 
If we let e -> 0, we firstly have that 



where tp s is the solution to i|8.20| with discontinuous initial data 2Ma(x - s); and 
secondly that for all t < r, 

u*(0,t) < ^(0,t), 
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where s is chosen so that tt(0, t) = <^ s (0, t). 

Now we turn to the general domain with Q = [x ,xi). Given i £ fl, set d := 

d 2 TT 

distfx, 90). If t < — ; , — -— r then we can repeat the same calculation on the small 
domain [x- d,x + d] to find the given result. □ 



Chapter 9 



Estimates for isotropic and 
anisotropic mean curvature flow 

9.1 A gradient estimate for mean curvature flow 

This section follows the style established in papers such as (T3) and (T4), in partic- 
ular the local gradient estimates of Section 2 of the latter paper. Ecker and Huisken 
consider the evolution of a hypersurface by mean curvature 

jF( P ,t)=H( P ,t), P eM, (9.1) 

where F : M n x [0, T] — >■ R n+1 is the immersion of the manifold M n at each time t and 
H is the mean curvature vector. 

M can be written as a graph when a fixed vector u e R n+1 can be found so that 
for a choice of unit normal v, 

(u,oj) > 

everywhere. Equivalently, (u,^)^ 1 is bounded above. The existence of an upper 
bound for this quantity and its analogue for anisotropic mean curvature flow is the 
subject of this chapter. 

Given the image F(p,t) of a point p e M, its coordinate vector is x(p,t). The height 
of M above the hyperplane defined by u is denoted by 

u = (x, (jS). 

The gradient function is given by 

v = {v,uj)~ 1 = a/1 + \Du\ 2 
We recall the evolution equations from [13]: 
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Lemma 9.1. lfM t satisfies (|9.1f then 

d 



A )|x| 2 = -2n, 

dt 



d 
d~t 



A)v = -\A\ 2 v -2v~ 1 \Vv\ 2 , 



where A /'s tf?e Laplace-Beltrami operator on M t and A = {h i:j } is the second funda- 
mental form. 

We derive a gradient estimate for periodic entire graphs, followed by an interior 
estimate. Estimates of this type have also been recently found by Colding and Mini- 
cozzi (T2) in the isotropic case using similar techniques, although without the explicit 
dependence on the height of the graph that the following estimates display. 

Theorem 9.2 (Estimate for periodic mean curvature flow). Let F be a smooth, en- 
tire solution to mean curvature flow (|9.1| which is a periodic graph over a hyperplane, 
in thatu(y, t) = u(y+L,t) for a fixed point L in the hyperplane, and has a height bound 
\u\ < M. Then 

forO <t < T', where c and T' depend on M. 
Proof: Define a new quantity 

Z := v — <p(u, t) 

where <p is a smooth positive function on [-M, M] x (0, X") chosen so that tp -> oo as 
t — > 0. This means that Z is strictly negative initially, regardless of the initial gradient. 
The evolution equation for <p is given by 

^ - A^ ip = ip t - ip uu \Vu\ 2 , (9.2) 

and we can use this with the identities from Lemma lPl to find that 

(jt - AJ Z = -\A\ 2 v - 2v~ 1 \Vv\ 2 + ^ u |Vn| 2 . 

Now, suppose (x,t) is the first point at which Z becomes non-negative. Since Z is 
periodic, this is an internal spatial maximum, and (spatial) first deriatives are zero: 

= VZ = Vv - Vip, 

so 

<p <p 
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and the evolution equation at this point is 



— - a] Z = -\A\ 2 v - 2^ |Vu| 2 -ipt + ^ u «|Vu| 2 . 



dt J tp 

A good choice for ip that will allow us to make the final terms negative is ip(u,t) 
l/$(u,t), where $ solves the heat equation & t = c<t>" for some c < 1. In this case 

tp> = -$-2$', 
(ff = 2$" 3 ($') 2 - $~ 2 $", 
ft = -^~ 2 ^t 

and the equation satisfied by ip is 

(ft = op — 2c . 

<P 

The final three terms of the evolution equation have become 

|Vn| 2 - ^ + ^ u |Vu| 2 = -2^ |V^| 2 + - p |Vu| 2 + 2^ |V^| 2 



§2 (c-|Vn 12 



In the last line we have used that, with respect to a local orthonormal frame on M u 

Vn = {ei,oj)eu while uj has unit length with 1 = \oo\ 2 = Ya=i I <e ; , c^>e £ | 2 + \{v,u))u\ 2 : it 
follows that 

\Vu\ 2 = f>,o,} 2 = 1 " ("M 2 = (l - ^) = (1 - $ 2 ) , 

the last equality holding only at a maximum point. 

If we let <J> be a fundamental solution of the heat equation 

t) = ^ exp f-ci ^ , (9.3) 

we can choose T and c depending only on M so that <t>" > and c - 1 + <3? 2 < for 

t < r. 

So, at the first interior point where Z = 0, Z t < and so Z < for t < T . □ 

Theorem 9.3 (Interior estimate for mean curvature flow). Let F be a smooth solu- 
tion to mean curvature flow (|9.1| which is a graph over a ball in the hyperplane B R (0). 
Then we have the interior estimate 

v < t«l 2 exp ( Cg(n + 2M)2 ) (R 2 - 2nt - |x| 2 + u 2 )" 1 
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forO<t<T', where q>l,c, and T' depend on M and R. 



Proof: We replace tp in our previous definition of Z by p/rj: 

<p(u,t) 



Z := v 



V 



where a smooth positive function rj is chosen so that Z < on the boundary of a ball 
of shrinking radius B ^ R i_ 2nt , and, as before, <p > is chosen so that Z < at the 
initial time. 

In particular, choose rj = R 2 - 2nt - |x| 2 + u 2 . The evolution equation for rj is given 

by 

|-a) ?? = -2|Vu| 2 , 



and we can use this, the identities from Lemma RTTl and the evolution equation for (p 
to find that 



^ - A J Z = -\A\ 2 v - 2v~ 1 \Vv\ 2 --(cp t - ip uu \Vu\ 2 ) 
at J rj ' 



- 2^ |V< - 2^Vu • V? ? + 2^|Vr ? | i . 

7]^ if 7]° 

Now, suppose (x,t) is an internal point of this domain at which Z first becomes 
non-negative. 

At an internal spatial maximum of Z, VZ = so 

|V,| 2 = v(^ 

= 4lV^| 2 -2^V W .Vr ?+ 4|Vr ? | 2 . 
r\ z rf rf 

Use this to replace the v^ 1 \Vv\ 2 term in the evolution equation, so that at this point 

d a\ 7 1/112 oVf&WTfl oW«U Y7 I V 2 IT-7 1 2 

- — A Z = — L4 v — 2— — £ Vw - 2 — r Vu • V77 H — j Vr? 
at J f \V V V 

--(<pt- ^ uu \Vu\ 2 ) - 2^ |Vu| 2 - 2^Vu ■ Vr, + 2^|Vr ? | 2 
= -\A\ 2 v --(cpt- <Puu\Vu\ 2 + 2 |Vu| 2 ^ 



2|Vd 2 -^ + 2^Vu- Vrj. 
rj z rf 



We can bound the Vu • Vr? term by a v~ l term, since (using a local orthonormal 
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frame {e^}) 



and 



Vu- Vt] = V(x,w) • V(-(x,x) + (x,w) 2 ) 

= (VjX, uj)ei • (— 2{x, Vjx)ej + 2(x,u)(VjX.,v)ej) 
= (ei,u)ei ■ (-2(x, ej)ej + 2(x, w)(ej, u)ej) 
= 9 v (ei,v) (-2(x,ej) + 2(x, u})(e jt u)) , 



l + \Du\ 2 
(ei,u) = D, t u 
(e*,x) = x» + 

so that (writing x for x - w, the position in the hyperplane) 

(Du,x)\Du\ 2 



Vu ■ Vr] = -2{Du,x) + 2 
= -2 
< 2 



l + |L>ii| 2 

(Du, a;} 
'1 + \Du\ 2 
x 



< 2 



y/l + \Duf 
VR 2 - 2nt 



v 

With this we estimate the term in the evolution equation — 

2^Vu • Vn < A^R. 
The evolution equation itself becomes 



|-a)*<-|^-2|V»|^ 



1 ( iv7 |2 , |v7 |2 ^ a r> I I 

<Pt-<Puu\vu\ + 2|Vu| 4R- 



This time, we choose 

<p{u,t) = <S>(u,t)~ q 
for some q > 1, where $ still satisfies the heat equation <t> t = c<t>". Then 

iff = -q®-*- 1 ®' , 
ip" = q (q + l)$-'?- 2 ($ / ) 2 - q<5>-i- l <&" , 
(ft = -qQ-t- 1 ^ 
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so that the equation satisfied by <p is 

1\ <// 2 



ip t = c<p - c 1 + 



The final term in the evolution equation is 

2 I I 

<Pt ~ Vuu\Vu\ 2 + 2 \Vu\ 2 ^ - ABV' 



/ ON l$'l / n |$'| \ 

= 9i5+T (" c + 1 " * 9 ) + 9 V ( (Q " " * " 4R ) ' 

where we have used \Vu\ 2 = 1 - = 1 - <3? 2|J . The first term above is positive if, 
as in the previous case, 3>" > and c — 1 + <t> 2g < 0. The second term is positive if g 
satisfies 

9 " 1+ |$'|(l-$ 2 )' 

Choose $ to be a fundamental solution of the heat equation, as in the previous 
proof — 

, . 1 / (u±2M) 2 
&{u,t) = —= exp —c- 



Now we can choose some c < 1 and T" small so that > and <J> < (1 - c) < 1 for 
t < T'. Here, T is dependent on M and c only. We can also find q dependent on T', 
M, c and R, satisfying 

1 + | $'|(1 -$2) " 1 + C M(1 - (1 - C) 2 ) - 9 ' 

At an internal maximum of Z t < and the result follows. □ 



9.2 Gradient estimates for anisotropic mean curvature 
flows 

A more general case of curvature flows is that of anisotropic mean curvature flow. 
This has been specifically studied by Almgren, Taylor and Wang (TJ, Gurtin and An- 
genent 0, and Andrews 0I31, among others. The anisotropic surface energy arises in 
applications from materials science, such as crystalline growth and phase changes; it 
also arises in Finsler geometry (Tl] (on a Finsler manifold, at each point only a normed 
space is defined, rather than an inner product space as on a Riemannian manifold). 

In this section, we use the framework and notation of |4). 

As before, we consider surfaces that can be written (either entirely or locally) as 
graphs, so that M t = {(x 1 , . . . ,x n ,u(x 1 , . . .,x n ,t))} = graph u(x,t). 

The equation for motion of the graph by anisotropic mean curvature is derived in 
0; in the present work, we set the homogeneous degree zero "mobility factor" m to 
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be identically 1, and so 

ut= F(Du)D^F\ DuUij , (9.5) 

where u„ = D 2 u(ei, ej), with respect to some basis for the tangent space {e 1 ,...,e n }. 

The function F : R n — >■ R is defined by F(p 1 , . . . ,p n ) ■— Fipift - 4>°), where 
{<f)°, (j) 1 , . . . , <j) n } is a basis for the cotangent space V*, with dual basis for V = R n+1 
given by {e , e±, . . . , e n }, and where F -. V* —> R is a positive convex function that is 
homogeneous of degree one, F(\u) = \F(lu) for A > 0. The unit ball of F 

F^il) := {w £ : F(u) = 1} 

must be strictly convex. We also require that F is at least C 3 . 

Differences between the isotropic and anisotropic cases 

The introduction of the unspecified anisotropic F into the flow has the effect of high- 
lighting the special nature of the isotropic case, when F(p i (f) i +po4>°) = {Ya^pI) 1 ^- 
One immediately notices that in the isotropic case, the term with third derivatives 
arising in the evolution equation is zero. In Q9.21| . this is the term 

DF\ z (0 fc ) D (FD 2 F) \ z (<T, 4>\ &)u mkUij . 
This absence of third derivatives is apparent in the third identity of LemmalO] — 

[j t -/^jv = -\A\ 2 v-2v~ 1 \Vv\ 2 ] 

the left-hand side involves second derivatives of the gradient so we might expect to see 
some derivatives of curvature in the right hand side — instead we see only curvature 
terms and first derivatives of the gradient function. 

The second difference is that there is no estimate of the form 

Vu ■ Vi] < -, 

v 

as in \§A) for the isotropic case. An equivalent estimate in the anisotropic context 
would be: given q = Yh=i f' nc ' c = c (l) so tnat 

F(p-<jP)FD 2 F\ p ^ {p,q)<c(q) 

for all p = YTi^Vi^- This is certainly true if we restrict p to the unit ball, F{p) = 1. If 
we replace p by sp (s is a scalar), the putative estimate would be 

F(sp - <f>°) FD 2 F\ sp ^ (sp, q) < c(q). 

Rewriting the left-hand side using homogeneity gives us 

sF(p- cjP/s) F£> 2 F|__^ 0/s (d>°,q) < c{q). 

As s increases, the left-hand side is converging to a constant defined on the unit ball, 

F(p) FD 2 F\. 0°, q), multiplied by s. Unless F(p) FD 2 F\. (<p°, q) is zero , the left-hand 
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side will not remain bounded by the right-hand side as s — ► oo. The estimate will not 
hold without further restrictions on F. 



Calculating with the homogeneous function F 

We make some observations about properties arising directly from the homogeneity 
and convexity of F, and introduce some notation. 

Let {<p°, . . . , 4> n } be a basis for the cotangent space V* dual to {e , . . . , e n }, the 
basis for the tangent space V. Both V* and V are copies of R n+1 . 

For p = picj) 1 (all repeated indices are summed from 1 to n unless indicated other- 
wise) we will write 

z := p — <p° . 

In general we will prefer to write all derivatives of F in a form that is homogeneous 
of degree zero, that is, as DF, FD 2 F, or F 2 D 3 F. This means that we can evaluate 
them on the unit ball, or scale as we wish — for example, we can use FD 2 F\ p _ (f>0/t 

instead of FD 2 F\ . 

Homogeneity also means that some derivatives in the radial direction disappear — 

for all uj £ V*, 

DF\ u (u>) = F{u) (9.6) 
FD 2 Fl (oj, •) = FD 2 Fl (;u) = (9.7) 
D(FD 2 F)\Jco,;-)=0. (9.8) 

The strict convexity of the unit ball of F means that for all u, r e V* on the unit ball, 
with r / ±uj, 

D 2 F\ (r,r)>0. 

As F is homogeneous, all the level sets of F are also strictly convex, so this holds for 
all non-zero uj. 

We denote by ~ the removal of a component in the direction of z from <j) k , k = 

0, ...,n, 

fc : = (j) k - c k z, 
where c k is such that 4> k is tangent to the unit ball of F, 

0=DF\ z ($ k ) 
= DF\ z {<f) k - c k z) 
= DF\ z (c/> k )-c k DF\ z (z) 
= DF\ z {<p k )-c k F{z). 

We have used i|9.6| in the last line. It follows that 

DFim 

F(z) ■ 
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In the next two lemmas, we show that the coefficients of the evolution operator 
satisfy a condition similar to the control on degeneracy that we required with condition 
(ED of Chapter^ 

Lemma 9.4. Let F : V* -> R be a C 2 , positive, homogeneous degree one function 
with a strictly convex unit ball F~ l {l) = {uj : F(uj) = 1}. 

Let 0°, . . . , <p n be a basis for V* . 

Then there exist positive constants A and P so that 

FD 2 F\ p _ (j>0 (p,p)>A 
for all p = ElLi^^ with F{p) > P. 



Proof: Write 

B(p) := FD 2 F\ p ^ (p,p) = FD 2 F\ p ^ (p,p), 

where 

DF \ m (p) 
p p F (p _ 00) W V ) 

is non-zero whenever p^O. As p is a non-zero tangent covector, the strict convexity 
of the unit ball means that B(p) > if p ^ 0. 
Fix p = picj) 1 on the unit ball of F, F(p) = 1. 

Consider B(sp). We want to show that we can find some P p and some strictly 
positive A p so that for all s > P p , B{sp) > A p . If this is not possible, then we can find a 
sequence s k — > oo with 

lim B(skp) = 0. 

fc— »oo 

Since F is C 2 , FD 2 F\ is continuous in z, and we have that 

Hm ™ 2 FL fcP „,o (-AP, s kP ) = Km ^| SfeP _,o (0°, /) 

" lim ^ 2 ^| p - o /sfe (0°,0 ) 



fc— >oo 

FD 2 F 



Clearly, 

lim p - 0°/s fc = p, 



A;— >oo 

SO 



lim S( Sfe p) = FD 2 F| = ^ 2 F| > 0, 



by the strict convexity of the unit ball, as at p, 



DF\ i 



F(p) 

is a non-zero tangent covector. The contradiction implies that we can indeed find such 

P p and A p . 
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We can find such P q and A q for every q = q^K Let 

A := jnf A q , P := sup P q . 

q--F(q) = l q:F(q) = l 

As we are optimizing over a compact space, A > and P < oo. 
The result follows directly, for given any p = pi4> 1 with F(p) > P, 

B(p) = B(F(p)p) > Ap > A, 

where p = p/F(p) is on the unit ball. □ 

We use this to show that the anisotropic mean curvature flow satisfies i|5.9| . the 
condition controlling the degeneracy of the parabolic operator in Chapter|5l 



Lemma 9.5. For all non-zero v = andp = p^, we can find positive constants P 
and Aq such that 

FD 2 F\ p _ (j)0 (v,v)>A (9.9) 



(v-p) 



whenever F(p) > P. Here, (v ■ p) 2 = Y^( v iPi) 2 an d \p\ 2 = p-p- 



Proof: Set B(p, v) = ^ FD 2 F\ p ^ (v, v). 

Since B(-,v) is invariant under v h-> sv, we need only to consider v in the unit ball. 

Suppose that p is in the unit ball. Since v is a non-zero tangent covector at p - <fP, 
FD 2 F\ p _ (f)0 (v, v) > by the strict convexity of the unit ball. By compactness, 

inf inf FD 2 F\ , (v, v) > ci > 0. 

P eF-i(i)„eF-i(i) 

Also, as neither p nor v are zero, 

M 4 



(vp)' 



> C 2 > 0, 



and so inf p veF -i^ B(p,v) > cic 2 > 0. 

Suppose, in order to obtain a contradiction, that there is a pair (v,p) in the unit ball 
for which there are no such constants A and P. That is, 

lim B(sp,v) = lim f FD 2 F\ , (v,v) = 0. 

s^oo V ' K ' ' s^oo( V -p) 2 \sp~<f> ^ ' I 

There are two possibilities here: v / p or v = p. In the first case, we must have 

= _lim FD 2 F\ sp _^ (v,v) 



♦oo 



lim FD 2 F\ 



= FD 2 F\ p (v,v). 
However, since v^p, FD 2 F\ (v, v) > which is a contradiction. 
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On the other hand, if v = p, then 

lira B(sp,p) = lim s 2 FD 2 F\ sp _ <j>0 (p,p) 



= }™o FD F \s P ^o( s P^P) 

> A 

by Lemma l9~4l which is again a contradiction. 

Therefore for every pair of covectors (v,p), there is a pair of positive constants A 
and P such that B(p, v) > A whenever F(p) > P. To get bounds for all (v,p) we take 
the infimum of the A and the supremum of the P. □ 

We will consider two different restrictions on F. The first is that third derivatives 
are small; the second is a symmetry in the distinguished direction 4>°. 

In order to define the first condition, consider the tensor 

Q(p, q, r) := F 2 (z) D 3 P\ z (p, q, r), (9.10) 

for p, q, r covectors tangent to the unit ball of F at z, so that DF\ z (p) = (and similarly 
for q and r). 

(This is the Cartan tensor of Bao, Chern and Shen [TD|. or the tensor Q of |4] 
restricted to the tangent space of the unit ball.) 

The smallness-of-third-derivatives condition is that Q satisfies 

Q(p, q, r) < d [F{zf D 2 F\ z (p,p) D 2 F\ z (q, q) D 2 F\ z (r, r)] 1/2 , (9.1 1 ) 

for all p, q, r tangent to the unit ball of F, where C\ is a positive constant dependent on 
n. 

The symmetry condition is that 

n 

F( P + = F(p - cjP) for all p = ^Pif. (9.12) 

i=l 

Lemma 9.6. If F satisfies (|9~T2| . then 

DF| p (0°) = O (9.13) 
D 2 F\ (^°,^) = 



I 

\p 



^F| p (^,^,^) = (9.14) 



^F| p (0 u ,0 u ,0 u ) = O, (9.15) 
for allp = Y%=i Pift and allj, k^O. 



Proof: This is a direct consequence of homogeneity. □ 

We can show that the symmetry condition (|9.12| can be used in a similar way to 
the smallness-of-third-derivatives condition (|9.1 1| . 

Lemma 9.7. Suppose the symmetry condition (|9.12| holds. Then a constant depen- 
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dent only on F bounds 



FD(FD 2 F) 



p-<P 



,o (p, q, q) 



(9.16) 



G(p,p)V*G(q,q) 

for all p = YTi^Vi^ andq = £™ =1 where G = FD 2 F\ _ 0O . Furthermore, for all 
e > we can find S t so that when F(p) > S e , 



FD(FD 2 F)\, (p,q,q) < eG(p,p) l ' 2 G(q,q). 



Proof: Consider p and q on the unit ball and set 

FD(FD 2 F)\ (p,q,q) 
v ,qZF-\\) G{p,py/ 2 G{q,q) 

p=Pi<t> 1 ,q=qi<j) 1 

When we project p and q onto the tangent plane atp-(p° they give non-zero tangent 
covectors p and q, 

p = p-( p (p-4>°), q = q-c q {p-4>°), 

where 

DF \ m(p) df \ m(q) 

F(p-(ffl) ' F(p-cJ) ) ' 

so the terms in the denominator of Q9.16| . G(p,p) and G(<?, g), are strictly positive, and 
hence bounded below when we take the supremum over p and q in the unit ball. 

Also, FD{FD 2 F) is a homogeneous degree zero tensor, and so bounded above 
on the unit ball. 

It follows that C is finite. 

The constant C is unchanged if we scale q h-> sq so we only need to consider the 
behaviour of (|9.16| as F(p) becomes large: that is, 

FD(FD 2 F)\ , (sp,q,q) 
s-»oo G(sp,spy/ 2 G(q,q) 



where G = FD 2 F\ 



sp—<jfi ' 



Firstly, consider the case that q is parallel to p. Let p be on the unit ball, and without 
loss of generality, let q = +p. 

We note that the G(q,q) = G(p,p) term in the denominator converges to zero, 



lim G(p,p) = lim FD 2 F\, ,(p,p) 



p—(jfl/s 



FD z F\ p (p,p) = 0, 



96 



so to deal with this we will multiply both top and bottom by s 2 : 

FDiFD 2 ^)]^ (sp,p,p) 



lim 



s^oo G(sp, sp) 1 / 2 G(p,p) 

s 2 FD(FD 2 F)\ gp ^ (<jP,p- C P(sp - <jP),p - C P(sp - /)) 



lim 



lim 

s— >oo 



s 2 G(sp, sp) 1 / 2 G(p,p) 
FD(FD 2 F)\ sp ^ <p0 (</>°,s(l - scP)p + scP(j) ,s(l - s C P)p + 8<PtfP) 



. (9.18) 



G(sp, sp) 3 / 2 

Now, the denominator is strictly positive, and by Lemma l9~4l bounded below — 

lim G(sp,sp) = lim FD 2 F\ sp _ A0 (sp,sp) > A > 0. 



The limiting value of the coefficient of cf>° in \9A8) is 



DF 



p-<t>°/ t 



(p) 



lim scP = lim s- - 

s^oo s^oo F{sp — (/)") 

DF \ AO/ (P) 

= lim _ 'p-* /« w 

s->oo F{p — (fP/s) 

DF\(p) 



F{p) 

= 1, 

using Q9.6| . 

The limiting value of the coefficient of p is 

lim s(l — scP) = lim s I 1 — s 



lim - 1 

r->0 r \ 



DF \ P ^/s (P) 

F(sp - 0°) 
^U^o(p)\ 



lim 



F(p - r0°) / 
1 fDF\jp) DF\ p ^ (pY 



•or 1 F(p) F(p - r</>°) y 
where we use that DF\ (p) = -F(p). The above term is a derivative, so we have 

d ( DF \ P +r4>° b>y 



lim s(l — sc p ) 



dr I F(p + r<; 



r=0 



D 2 Fl (0°,p) -D-F^L (p) DF\ (4>°) 



F(p) 



Hp) 2 



0. 



where the first term of the second last line is zero due to i!9.7t while the second term 
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is zero as consequence <9.13| of the symmetry condition. 
The equation (|9.18| is then 

FD(FD 2 F)\ s , (</>°, s(l - sc p )p + sc p </>°, s(l - sc p )p + sc p (/) ) 
lira ^— - 



G(sp, sp) 1 / 2 G(sp, sp) 



FD(FD 2 F) 




FD 2 F 


p (0O ) O )3 /2 



= o, 

where symmetry has been used again in the form of (|9.15| . 

Now consider the case that q is not parallel to p. In this case, the denominator of 
i|9.17| is non-zero: 

Urn G(sp, sp)^ 2 G(q, q) = Km FD 2 F\ p _^ 0/s 0°) F£> 2 F| p _^ 0/s ( g , g ) 

= FD 2 F\ p (<l> ,<f> ) FD 2 F\ p (q,q) 
>0, 

since lim s _*oo g is a non-zero tangent vector at p, 

DF| DF| (<?) 

lim ~ 9 = lim ~ 9 - F( S p~- (sp " * 0) = 9 " ^r p - 



■>oo s— >oo 



This cannot be zero as g is not parallel to p. 
Then 

s^oo G{sp,sp) l l 2 G{q,q) 

-1 



= lim ( G(sp, sp) 1 ^ 2 G(q, g)J lim FD{FD 2 F)\ p , 0/s (^,q,q) 

= o, 

where the second limit is zero by i|9.14fr . since lim^oo q has no component in the 
direction of 4>°. 

We have shown that for a fixed p on the unit ball, the quantity \9.~\6\ is bounded 
above (by C), and that as p is scaled outwards this decreases to zero, so for fixed p 
the quantity is bounded above. 

By compactness of the unit ball, it is bounded for all p. □ 

Estimates for periodic, anisotropic mean curvature flows 

Let u : R n x [0, T] -► R be a H 2 , bounded 

\u(x,t)\ < M, 

periodic 



u(x + L,t) = u(x, t) for some lattice L, 
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solution to the anisotropic curvature flow equation i|9.5| . where F is a positive convex 
function, homogeneous degree one, with a strictly convex unit ball. 

Remark: We have one estimate in the case that F satisfies the smallness-of-third- 
derivatives condition (Theorem I9.8| and another in the case that F satisfies the sym- 
metry condition (Theorem I9.9|. However, in Theorem l5?TI and Corollary l5^2l we found 
an estimate for periodic anisotropic flows without the need to impose either of these 
conditions. In that case, the fact that we were estimating on the difference quotient, 
rather than the first derivative itself, avoided the need to take derivatives of the flow 
coefficients. Strictly speaking, Theorems MM and |9JD are redundant, but they are a 
good introduction for the interior estimate of Theorem r9.12l 

Theorem 9.8. If the tensor Q given by i|9.10| satisfies (|9.1 1| with 

Ci 2 <^=, (9.19) 

then 

F(^)< m ax{^exp(^H-2Mfy p | 

for < t < T , where T depends on M, A and P (both given by Lemma \5Hj> . and 

1 < q = (1 - Cfv^/4)- 1 . 

Theorem 9.9. If F satisfies the symmetry condition i|9.12f . then 

F(Du) < maxjtexp , P, S (2/n)1/3 } 

for <t <T', where T' depends on M, A and P (both given by Lemma WM . and S is 
given by Lemma \9l\ 

Proof of Theorem l9.8t As in the previous sections, define the quantity 

Z := F(Du) - <p(u,t) 

where P is given by Lemma l9~4l and where <p is a positive function chosen so that 
ip oo as t -> 0. Suppose that we are at the first point where Z is no longer negative, 
and let us assume that at this point, tp > P. This point will be a spatial maximum of Z, 
due to the periodicity of Z. 

Then at this point, F(Du) = (p and the first derivative condition is 

= D k Z = DF \ z ((j) m )u mk - ip u u k , 

where z = Du- <jP. That is, for all vectors v e span{ei, . . . , e n }, 

D 2 u(DF\ z (<i) m )e m ,v) =<p u Du(v). (9.20) 

Using Q9.7| . we can rewrite the evolution equation for u in terms of the purely 
tangential directions 4>\ 

u t = FD 2 F\ z ^\^)D 1 u[e i ,e j ) = FD 2 F\ z {^ ^)D 2 u(e i ,e j ). 
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We make use of this in finding an evolution equation for F 



BF 

— = DF\ z (<t> k )u kt 



DF\ 



DF\ 



FD 2 F\U\4P)u ij 



D [FD F) \ z (D k z, 4> l ,^)uij + FD z F\ z i/ W^ i'',, 

+ FD 2 F\ z (4>\D k ^)u i3 + FD 2 F\ z ($\ft)u ijk 
+ FD 2 F\ z (jj)DijF 

- FD 2 F\ z {<j>\ </>>) [D 2 F\ z (</> m , 4> l )u mi uij + DF\ z (<p m )u mij 

= DF\ z (<p k ) [ D {FD 2 F)\ z (D k z, ft, ft) Ulj + FD 2 F\ z {D k ft, ft>) Uij 

+ FD 2 F\ z (ft,D k ft)u i3 
+ FD 2 F\ z {ft, ^ i )D lj F - FD 2 F\ z (ft,c^) D 2 F\ z (<j> m , ft)u miUlj , 

where in the third step we have added and subtracted second derivatives of F. 

Derivatives of z are D k z = u mk (p m . We this to simplify those terms with derivatives 

of ft— 

D 2 F\ z (D k ft,ft) = D 2 F\ z (D k (-c^z),ft) 

= D 2 F\ z (-D k (c^)z - c^u mk( t> m ,ft), 

and remembering that D 2 F\ (z, •) = 0, this is 



-<f n 2 F\ z (u mk <r^) 

DF\ (ft) „_. 

- D 2 F\ z {u mk r^: 



The evolution equation is now 



dF 
~dt 



DF\ (<f) k ) D (FD 2 F) 



DF\ (cf> k ) DF\ (ft) D 2 F\ 7 (<T, ft) + DF\ (ft) D 2 F\(<p m , ft) 



'U"mk u ij 



FD 2 F\ z (f,^>)D lj F - FD*F\ z (<?,</>>) D 2 F\ z (^ m ,ft)u mi u lr (9.21) 



When we are at a critical point of Z, we can use the first derivative condition i|9.20> 
to simplify further. The first term of (|9.21| becomes 

D 2 u(DF\ z (<f> k )e k , e m ) D (FD 2 F) \ z (4> m , ft, ft) Uij 

= tp u Du(e m ) D (FD 2 F) 



^,<P\^) Uij 



Vu D(FD 2 F)\(Du,ft,ft)u. 
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while the second becomes 

- DF\ z {ft) DF\ z {ft) D 2 F\ z {<j> m , ft)u mk u i3 

= - D 2 F\ z {<j> m , ft)D 2 u (DF\ z {4> k )e k , e m ) D 2 u (DF\ z (ft) e h ej ) 

= -ip 2 u D 2 F\ z {Du{e m )0 m ,Du{ ej )ft) 
= -iflD 2 F\ z (Du,Du), 

as does the third. 

The evolution equation for F, at the local maximum of Z, is now 

^L = ^FD {FD 2 F) \ z {Du, ft, ft) Uij - 2^ FD 2 F\ z {Du, Du) 

+ FD 2 F\ z {ft,^)D ij F-- FD 2 F\(ft,ft) FD 2 F\U m ,ft)u miUl] , 
' z tp z 

where we have multiplied some terms through by 1 = F/p (since we assume that 
Z = here) in order that derivatives of F appear as homogeneous degree zero terms. 

Derivatives of p are given by 

Dp = p u Du 

Dijtp = ip uu UiUj + tp u Uij 

dp 

— = p u u t + tp t 
at 

for i, j / 0, so an evolution equation for tp is 

^ = ip u U t + p t + FD 2 F\ z ft) {Dijtp - PuuUiUj - tpuUij) 

= pt + FD 2 F\ z {4>\ ft) {Dijp - p uu UiUj) , 
and the entire evolution equation for Z, at a local maximum, is 

^ = FD 2 F\ z {ft, ft)D ii Z + ^FD {FD 2 F) \ g {Du, ft, ft) Uij 

_ FD 2 F\ z {Du, Du) - - FD 2 F\ z {ft, ft) FD 2 F\ z (<T, ft)u miUlj 
p ,z p ,z ,z 

-p t + FD 2 F\ z {Du,Du)p uu . 

Notice that all the covectors <fr\ Du appear in places where they may be replaced 
by ft, Du respectively (using i|9.7| and (|9.8|). That is, we are working exclusively on 
the tangent space to the unit ball. 

Restricted to the tangent space, D 2 F is positive definite so j/ve can define a Rie- 
mannian metric on the tangent space G af3 := FD 2 F\ z {4> a ,4> 13 ), a,/3 / . We 
can choose the basis {ft, . . . , 4> n } so that G is the identity at our maximum point, 
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Qo.fi _ j^g evo | u tion equation for Z is now 

— = G ij DijZ + ^FD (FD 2 F) \ g (Du, fi^Uij - 2^G(Du, D 



dt J <P K ,T,T " J tp 

G^G ml u miUlj -(pt + G(Du, Du)<p uu . (9.22) 



Recall the Cauchy-Schwarz inequality for positive matrices: If A is a positive semi- 
definite n x n matrix, then for v,w e W 1 , 

< ((2e) 1/2 v - (2e)- 1/2 w^ A ((2e) 1/2 v - (2e)^ 1/2 

and so 

v Aw < ev Av H w Aw. 

4e 

If A is positive definite, then we can replace w by A~ l w to find that 

v w < ev Av H u; A~ w. 

~ 4e 



As we assume that u is smooth, G is positive definite and we can use the above 
inequality to estimate the second term of \9.22) : 

— FD (FD 2 F)\ z (Du,fi,ft) Uij 



fu 



DF\ (Du)FD' 2 F\ (c/)\ (/>>) + F 2 D 3 F\ (Du,4>\ 



u 



i.) 



< e^G(Du, Du) + -^-G^Q^UijQ^Uku (9.23) 

where the first term of the second line is zero because Du is tangent to the unit ball, 
so DF\ z (Du) = 0. In the last line, we have used the notation for the inverse G a/3 = 
(G- 1 )^ 3 . We will choose e < 1 later. 



We can use \9.~\1\ . the smallness-of-third-derivatives condition, to estimate the 
second term in this inequality: 

l-G^Q^u^uu = -±-Q (G a ^ a , Uij f,^ Q (jPtUuPrf) 

cl 



^ T~A G(G a p<r,G^)G(u i ^\u m3 ^)G(^>^) 

_ _ _ \l/2 

x G(^, <f)G(u kl <p\ UpK/fMt 1 ^ 1 ) 
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= ^ (G aP G^G lP u ijUmj G^G^G^G^u kl u pl G li y 
= (GppUijUijG^G^uum^G 11 ) 1 ' 2 
= ^V^{G^G kl u ikUji y 



Now we can estimate \9.22) from above — 

r!7 C 2 2 

— < G ij DijZ + ^-^G ij G kl u ik u n + e^G(Du, Du) 
at Aeip tp 



2^G(Du, Du) - -G ij G ml u mi u l:i -tp t + tp uu G(Du, Du) 



tp tp 



G ij DijZ + 1 (^V^ - l) ('"•'(^ »//.";/ 



2 

+ ^ (e - 2) G(Du, Du) -tp t + tp uu G(Du, Du). 

The second term is zero if we choose 

e = d 2 V^/A < 1; 

the inequality here is a consequence of <9T9j. 

As in the proof of Theorem l9.31 choose tp = <t> 9 for some q > 1, with <t> given by 
<|9.3p> . This satisfies the heat equation <t> t = c$". We will choose c = A, where A is 
given by Lemma i|9.4Jl Substituting and its derivatives for tp and its derivatives (see 
page EH} we find that 



2 

^- (e - 2) G(Du, Du) -tp t + tp uu G(Du, Du) 



$-9-^'^( e - 2)G(Du, Du) 



+ 



]{q + \)<5>- q - 2 & 2 - q^- 1 ^" G(Du, Du) + Aq$~ q ~ l §" 
-2^/2 



g $-9-2$'^ [ g ( e - 1) + 1] Du) + g^" 9 " 1 ^' [A - G(Du, Du)} . 



Now, the first term is zero if we choose 

1 

Q 



1-e ~ l-Ci 2 Vn/4" 
As we assumed at the beginning that F{Du) > P, Lemma|SL4|ensures that 

G(Du,Du) > A. 

Consequently, if we choose T > small, is positive and then the final term will be 
negative. 

On the other hand, if we consider the possibility that tp < P at this local maximum, 



Chapter 9: Estimates for isotropic and anisotropic mean curvature flows 



103 



we could replace ip by sup{</?, P} in the definition of Z. In that case, the first maximum 
of Z occurs at a point where the barrier is flat, and so the first variation is 

= D k Z = DF\ z (ft)u mk , 

and the evolution equation for Z at the local maximum is 



dZ - o - 1 
— = FD F\ 
dt 1 

< 0. 



\ft>)D ij Z- D 2 F\(ft,ft>) FD 2 F\ z (ft n ,ft)u miUlj 



Since Z t < at the first point where Z = 0, Z < for all t < T and the conclusion 
follows. □ 

Proof of Theorem I9.9t We begin by assuming that at a local maximum point 
F(Du) > max{P, SV 2 / n )i/2}, and follow the proof of Theorem l9l3l up to equation i|9.22| . 
the evolution equation for Z at a local maximum: 



dZ 

~d~t 



G'WijZ + — FD (FD' 2 F)\ z (Du, ft, ft>) Uij - 2^G(Du, Du 



r 



1 



—G %0 G u m iUij -(pt + G(Du, Du)p v 



This time we do not choose coordinates to make G the identity 

We use Cauchy-Schwarz (with e = 1/2) to estimate the second term — 



FD (FD 2 F)\ (Du, ft, (jP)ui 



<Pu 



1 



< — G(Du, Du) + —G(Du, Du)' 1 FD (FD F) \ g (Du, ft, ft)uij 

Choose the basis {ft, . . . , 4> n } at this point so that D 2 u is diagonal. Also, let I 
be the n x n diagonal matrix with +1 or -1 as its diagonal entries, chosen so that 
F l uu = \uu\. As I 2 is the identity matrix, = F k \u k j\. 

In these coordinates, 



1 

2p 



G(Du,Du)~ L FD (FD F) (Du,ft ,ft)u { 



< 



< 



2<pG(Du,Du) 

1 1 

2^G(Du,Du) 

1 1 

2^G(Du,Du) 

n 



£ FD (FD 2 F) \ z (Du, ft, ft)f l \ Ul 



li=l 

n 



Y,\FD(FD 2 F)\ z (Du,ft,, 



1 

wp 



i=l 



£ 

i=l 
2 



U: 



1/2 



G(Du,Du) 1/2 G ii 



i=l 
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where the term (2/n) 1 / 2 comes from the estimation of FD (FD 2 F) | (Du, <fr\ in 

Lemma l9~7l under the assumption that F(Du) > 5 , (2/n) i /2 . 

This term is dominated by the fourth term of the evolution equation, since (in the 
same coordinates) 



-G 13 



ml 



1 



ml rma\ 



u miUlj = —G %0 G I 
<P 



1 

<P 



G ij 



jam Qml jl/3 



<p 



-G^G af3 \u ai \\ Ul3j \ 



1 1 

> — 

(p n 



'■j 

where we use the trace inequality (traced) 2 < ntrace(^4 2 ) 

2 



1 1 

if n 



in 



What is left of the evolution equation is 

^ < G lj DijZ - ^^G(Du, Du) -<f t + G{Du, Du)f uu . 

This is negative at a local maximum if we make the same choice of barrier as 
before — f = <&~ q for q = 2 with $ given by (|9.3| . 

If our assumption that F(Du) > max{P, 5 , (2//n )i/ 2 } does not hold, then we can 
replace if by max{P, 5 , ( 2//n) i/ 2 , if}. At the local maximum, Z t < and so the conclusion 
follows. □ 

Remark: In the last theorem, we have chosen q = 2 somewhat arbitrarily; in fact q 
needs only to be strictly greater than 1, since we can set q = (1 - ne 2 /4) _1 , for e given 
by Lemma I9T71 However, a smaller e may entail a larger S € , so the optimal choice 
would depend on the exact form of F. 



Interior estimate for anisotropic mean curvature flow 

We begin by showing that when we have the symmetry condition, an estimate analo- 
gous to (|9.4| in the isotropic case is possible. 



Lemma 9.10. Suppose that F satisfies the symmetry condition Q9.12| . Then there 
exists a constant C 2 depending only on F such that 



FD 2 F\^(p,q)<C 2 



F{p-<p°) 



for allp = p^ and q = qi 
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Proof: This estimate is unchanged under q i-> sq, so we need only to show that this 
holds for q on the unit ball. Let q = q^ 1 be a fixed point on the unit ball. 

By compactness, the estimate holds for all p on the unit ball. 

Let p be a fixed point on the unit ball and consider 

lim F(sp - 0°) FD 2 F\ (sp, q) 

= lim F(p-s0°)-FD 2 F\ (<jP,q) 
= YimF(p - s4>°)- s [FD 2 F\ p ^ (0°,q)- FD 2 F\ p (4>°,q) 

where we have added zero in the form of FD 2 F\ p (0°, q) 



-D(FD 2 F)\ p 



= F(p) l 

So, either the supremum of F(sp - cjp) FD 2 F\ (sp, q) over s e [0, oo) is the 

I sp <p 

limit above, or else it is attained at some finite s. In either case, 

C p (q) = snpF(sp - 4>°) FD 2 F\ sp _^ (sp, q) 

is finite, and we can set 

C2 = sup sup C p (q) 

{q:F(q) = l} {p:F(p) = l} 

to complete the lemma. □ 

The next lemma shows that the trace of FD 2 F is bounded below. 

Lemma 9.1 1 . Let {4>°, </>\ • • • , </> n } be a basis for V*, where n > 1. Then there is a 
constant k > so that for all p = Y^ =1 Pi^, 



1=1 

where G = FD 2 F. 

Proof: By compactness and strict convexity of the unit ball, 

sup Gl^p (<?,</><)= sup G\ p _^ o (ft,ft)>0, 

{p:F(p)<l} {p:F(p)<l} 

since ft is a non-zero tangent covector. 
Now consider 



t— >oo 

i=l 



where, in the limit, ft = 0* - D ^ ( ^ > V 
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At most one of the ft may be parallel to p — suppose it is ft, in which case 
lim^oo G\ t (ft, ft) = 0. For the remaining (n - 1) basis covectors, ft are non- 
zero tangent covectors (to the unit ball at p) and so G (ft, ft) is again bounded below 
for % / 1. 

It follows that 

n 

limY,G\ tp _0 Q (ft,ft)> G\ p (ft,ft)>0. 

i=l i=2,n 

If we take the infimum of all such lower bounds, over all p in the unit ball, then the 
conclusion follows. □ 

Let u : B R (0) x [0, T]^Rbeaff 2 , bounded 

\u(x,t)\ < M 

solution on the ball of radius R to the anisotropic curvature flow equation (|9.5| . where 
F is a positive, convex homogeneous degree one function, with a strictly convex unit 
ball. 

Theorem 9.12 (Interior estimate for anisotropic mean curvature flow). If F satis- 
fies both the smallness of third derivatives condition <9.1 1| with some constant 

C\ < 2/^i, 

and the symmetry condition \9A2) . then 

F(Du) < W 2 exp ^ q j\u\-2M? \ {r2 _ ^ _ |x|2) _, 

for <t< T'. 

Here, A is given by Lemma W^ and depends on F; k is given by Lemma\9jl\and 
depends on F; and T' > 0, q > 1, and r > 1 depend on M, A, P (which is also given 
by Lemma \9l\) and k. 



Proof: We introduce the localising term r/ into our definition of Z, now restricted to 
the shrinking ball: 

Z := F(Du) - - 
V 

for (x, t) g B j R i_ 2kt x [0, T], where k is the constant given by Lemma lPTI <p = ip(u, t) 
is a smooth strictly positive function chosen so that Z < at the initial time, and -q is a 
smooth positive function chosen so that r;^0on the boundary of the shrinking ball. 

Assume that at the first interior point where Z = 0, F(Du) > P. 

Then F(Du) = (p/r/ and as this is a spatial maximum (since the choice of -q ensures 
that there are no boundary maxima) we have a first derivative condition 



= D k Z = DF \ z (c/> m )umk ~ D k (if/r]) 



(9.24) 
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An evolution equation for <p/r) is: 



dt \r) ) 77 



FD 2 F\(4>\^) 



1 



+ 2^D i r } D^r ] - ^D ijV 



+ 2^G(Du, Drj) - 2^G(D V , Drj), 
where in the last line we have used the notation G ij = FD 2 F\ z (</>\ <^'). 



We can incorporate the first derivative condition \9.24\ into i|9.21| . the evolution 
equation for F: 

dF - 

— = FD 2 F\ z {if, (pP)DijF + DF\ z (</> fe ) D (FD 2 F) \ z (<j> m , f, &)u mk u i3 



- DF\ z (cf> k ) [DF\ z D 2 F\ z (<f> m ,ft) + DF\ z (^) D 2 F\ z ^) 
G ij DijF + D m D {FD 2 F) \ z (<T, ^')^ 



2^-G{D(<p/r ] ),D(<p/r ] )) - ^G ml u miUlj 



CPWijF + ^FD (FD 2 F) \ g [d^/t,), $ l , ^) u tj 



2 2 

^G(Du, Du) - 2^G(Du, Dr,) + ^jG (Drj, Drj) 
Tf rf r/ 4 

- ^G^UMj. 



Putting the last two steps together gives an evolution equation for Z at a local 
maximum: 



dZ _dF d (p 
dt dt dt \r) 

= G ij DijZ + ^FD {FD 2 F)\ z (z%/ ?? ) 



^j Uij -^G^G ml u miUlj 



ft ~ G(Du, Du)p> uu + 2-^G(Du, Du) 



(p ( d 



+ ^\Jt~ Gl3Di ^ rl+2 ^ G{DU ' Drl) 
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The second term here may be split up into a part with Dip and a part with Drf. 

^ FD (FD 2 F) | ( ^Du - ^D V , %\ ft) Uij 

= ^FD {FD 2 F)\ z (Du, 4>\ ~ X -FD {FD 2 F)\ z [dtj, u iy 

These may be individually estimated using the Cauchy-Schwarz inequality and the 
smallness-of-third-derivatives condition, as described in the proof of Theorem l9.8l on 
paae lToTI — 



FD(FD 2 F)\ z [Du 1 cp i ,(p j )u 



< W — G(Du, Du) + J-l(%yfti (G^G ml u miUlj 
1 FD (FD 2 F) I (dij, 4>\ ft) u 



< u 2 ^G(D v , Drj) + -L^C?^ (G^G kl u ih uA 
rr 4uo id V / 



Au 2 if 

for some < ux,u 2 < 1. 

We choose the localising term to be rj ■— ff for some r > 1, fj = R 2 - 2kt - \x\ 2 , 
and k given in Lemma|9jT| Then 

DiTj = rff~ l Dif) 
DijT) = rff~ l Dijf] + r(r - l)ff~ 2 DifjDjfj, 

and the second-last term of the evolution equation is 

^ (Jit ~ GijDi l) v = ^ r ^ _1 ^ 2k ~ 2traceG - ( r - l )v~ X G(Dfj,Df))] 
< ^rfj r - 2 (l - r)G(Dfj,Dfj). 

As F satisfies the symmetry condition i|9.12| . we may use Lemma lOOl to estimate 
the final term of the evolution equation: 

2^G(Du,Dr,) = 2^ FD 2 F\ Du _^ (Du,D V ) 
<2 Vu C 2 F(D V ) 



rj 2 F(Du - 0°) 
2C 2 F(D V )^. 



The evolution equation can now be estimated from above — 



^ < G ij DijZ + 1 (^-C 2 ^i + ^-C\^i - l) G^G ml u miUlj 
dt cp \4/xi 4u 2 J 
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2 

(ft - G(Du,Du)ip uu + (2 - iii)—G{Du,Du) - 2C 2 F(Dr))^ 

<P <P 

+ ^rff~ 2 {l-r + r^ 2 )G{Dfi 1 Df l ). (9.25) 



Since Cf^/n/A < 1/2, we can choose m < 1 and ^ 2 < 1 such that 

C%y/n (1 1 



+ - I < 1. 

With such choices, the second term of the evolution inequality <9.25| will be neg- 
ative. We can also set r = (1 - ^ 2 y l > 1, so the coefficient of f)~ 1 G(Dfj,Dfj) is 
zero. 

As in the previous cases we can set ip = where $ is given by \9.3) with c = A, 
where A is given by Lemma l9~4l 

The bracketted part of the second line of the evolution equation is then 

2 

tpt - G(Du, Du)(p uu + (2 - nx)—G{Du, Du) - 2C 2 F(Dr])^ 



9 



-q^-i- 1 ($ t _ G{Du, Du)<S>") 



+ G(Du,Du)q<f>' 2 $- q - 2 



-l + (l-/ii)?-2 



C 2 F(Dr]) $ q+1 
G(Du,Du) \¥\ 



(9.26) 



If we choose V small enough that <t>" > 0, then the term <t> t - G(Du, Du)®" 
(A - G{Du, Du)) is negative. 

As r > 1, F(Drj) = rfj r ~ 1 F(Dfj) < C s rR 2r ~ 1 , where C 3 > depends only on F. 
With this choice of <E>, 

2 _C 2 F{Drj)_^_ < 2 C 2 C 3 ri? 2r - 1 



G{Du,Du) |$'| 



< 4 



< 4- 



A 2 M 
,C 2 C 3 rR? r -H 



A 2 M 



if we choose T' small enough that $ < 1. 

In order to ensure that the last part of (|9.26fr is positive, we choose q so that 



> 



i - m 



1 + 4 



C 2 C 3 rR 



2r -\rpl 



A 2 M 



So, at such maxima, Z t < 0. 

At local maxima where F(Du) < P, then in the definition of Z we replace ip/r] by 
max{ip/r],P}, in which case the barrier is flat at the local maxima, and we again find 
that z t < 0. 

In either case, the maximum principle ensures that Z is never greater than zero 
and the conclusion follows. □ 



Appendix A 



Function spaces and regularity 
estimates for parabolic equations 



Here, we define relevant function spaces, and survey some regularity results used in 
the existence theorems of Chapters |U |6] and H This treatment follows the books of 
Krylov |23] and Lieberman |25) . 



A.1 Function spaces 

On the space of continuous functions u : Q — ► R, we have the supremum norm 

\u\0;U ■= SUp \u(x)\. 

Define the Holder semi-norm with exponent a e (0, 1] by 

\u{x) - u(y)\ 
[ujojn := sup — — . 

x,yeQ \ x ~ D\ 

For an integer k > 0, we define the Holder (k + a)-norm by 

\ u \k+a;Q ■ = 

|/3|<* \P\=k 

where /3 is a multi-index — an n-tuple of non-negative integers with \/3\ = Y^Pu and 

where D^n := — 1x „ — — -g-. 

The Banach space associated with this norm is C k+a (Q) = {u : \u\ k+a < oo}. 



The parabolic Holder spaces H a 

With parabolic equations, it is useful to weight the space and time variables differently- 
that is, two space derivatives to one time derivative. Following Lieberman, we will 
denote parabolic Holder spaces by H rather than C. For points z = (x,t) in a domain 
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Q c W 1 x R, define the parabolic Holder semi-norm by 

\u{zi) - u{z 2 )\ 

Na,a/2;Q : = S UP ; ; J— , 

2i,2 2 eQ (|xi - x 2 \ + \h - t2\ 2 r 
where a e (0, 1]. The parabolic norm is 

\ u \a,a/2;Q '■= \ u \o;Q + Ma,f ;Q- 

Higher spatial derivatives and derivatives in time are bounded by k + a norms, 
where k > is an integer: 

\u\k+ a ,a/2;Q = SUP \ D £ D t U \ + S [■ D * i> t U ]a,a/2- 

|/3|+2j<fc |/3|+2j=fc 

The Banach space associated with the | • \ k + a , a /2 norm is 

H k+a (Q) = {u : \u\ k+ata/2 < oo}. 

When the region Q is a cylinder, in the sense that Q = Q x [0, T], the parabolic 
boundary is given by 

v (n x [o, r]) = n x {o} uffix [o, r]. 

On the boundary, we can define parabolic norms H k+a (V) exactly as above. 

A.2 Regularity estimates 

In the following, P is a quasilinear parabolic operator 

Pu = —u t + a tJ (Du, x, t)DijU + b(Du, u, x, t) 

with positive constants Xk and A K such that 

a ij (x,p)^i> A^|e| 2 for^EM n 
la^'Cx.p)! < A K , 

whenever \p\ < K. We work on a domain O x [0, T] for some smoothly bounded 

n c M n . 

Here Q r is the intersection of the region and a cylinder: 

Qr{x , t ) := {(x, t) e n x [0, T] : \x — x \ < R,t - R 2 <t < t }. 

We begin with an oscillation estimate for the gradient of a solution for a Neumann 
problem near a flat boundary: 

Theorem A.1 . Suppose that a ij = a ij (Du, x),b = b(x), and that inside Q R , the bound- 
ary ofVt is x n = 0. Let u g C 2,1 n Hi(Q R ) be a solution of Pu = when x n > and 
D n u = 0, when x n = 0. 
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Suppose there are positive constants b and A such that 

\a ij p \ < A , \b{x)\ < b 

forall(x,p) e Q R x R n with \ P \ < K. 

lf\Du\ < K, then there are positive constants 9 and a determined only by K, n, A, 
A and A such that 



implies 



osc a %3 (-,p) < a for all \p\ < K 



osc Du < C(K, n,X, A, A ) ( [osc Du + b R\. 
nnQ r ~ \RJ nnQ R 



Similarly, we can find an estimate near the boundary for problems with Dirichlet 
boundary conditions: 

Theorem A.2. Suppose that a ij and b are uniformly continuous, that a ij is differen- 
tiate with respect to (p,q,x), and where, if \q\ + \p\ < K, we can find a positive 
constant such that 

K[\a» x \ + \a i3 q \\p\} + \b\<HK. (A.1) 

Ifu e C 2 ' 1 n Hi(Cl x [o, T]) satisfies Pu = 0onQx [o, T] andu = 0ondQx (0, T), 
then there are positive constants C and a depending onn,X K ,A K and ^kR/K such 
that for any (x , t ) efflx (0,T) 



osc Du < C 

Qr{x ,t ) 



— j K + +u K r 



forO < r < R <t 



1/2 



On the interior of the domain, one can also find a Holder bound for the gradient: 

Theorem A.3. Suppose thata ij and a are continuous; a ij is differentiate with respect 
to (p, q, x, t); and where if \q\ + \p\ < K, we can find a positive constant u K satisfying 

If u g C 2 ' 1 satisfies Pu = and \Du\ +\u\ < K in ft, x [0, T], then there is a 
positive constant a determined by n, Xk, A^,sup \a lj p \K, and u K R/K such that for 
interior sets x [h,t 2 ] cc fix [0,T] we have 

[ D u]a;U'x[t u t 2 ] - c ( n ' K i A^, A^,pt^,diam(r2 x [0,T]))d~ a , 

where d = dist(0' x [t u t 2 ],VVL x [0,T]). 

If Q r is a cylinder in the interior of the domain, we also have 

osc Du < C(n, K, Xk, ^k) (^) I osc Du + ukR 

Q r (x ,t ) \R/ \Q R (x ,t ) 

as long as < r < R < d(x ,t ), where d(x ,t ) is the distance from (x ,t ) to the 
parabolic boundary T>(tt x [0,t ]). 



Appendix A: Function spaces and regularity estimates 
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The following interior estimate is a H 2+Ct bound for u when the coefficients of P 
are smooth: 

Theorem A.4. Suppose that a ij = a ij (x, t), and b(Du, u, x, t) = b*(x, + c(x, 
a/id ffraf there is a constant K such that 

kM Q ,a/2 < K, 

for a e (0,1). 

Then for any R > fftere is a constant C dependent on R, X, a, K and n such that 

ifu e H 2+a (Q 3R ), then 

\ u \2+ a ,l+a/2;Q R < C {\Pu\ a ^ a/2 .Q 2R + Mo;Q 2H ) • 

In a similar vein, there are higher regularity estimates on the interior of a domain: 

Theorem A.5. If we have \D a a\ a ^ a/2 . Q2R , \D a b\^ a/2 . Q2R , \D a c\^ a/2 . Q2R < K for any 
|o| < k, and ifu e H 2+ct (Q 2R ) and D a {Pu) e H a (Q 2R ), then D a u e H a (Q R ) for 
\a\ < k, and there is a constant C = C(R, K, a, k, d, A) si/cft fhaf 

+a,l+a/2;Ofi ^ E l- D ^l2+a,l+Q ! /2;Q 2i? + CMo;Q 2jr 
|/3|<fe \p\<k 
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